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FOREWORD 


This report was prepared for the NASA Lyndon B. Johnson Space 
Center, Houston, Texas, by the Rocketdyne Division of Rockwell 
International. The study was conducted under contract NAS9-9S28, 
Mr. M. F. Lausten, Project Manager. Phases III and IV of a four- 
phase effort are reported. 


ABSTRACT 

This report contains the results (Phases III and IV) of a cold- 
flow and hot- fire experimental study of the mixing and atomization 
characteristics of injector elements incorporating noncircular 
orifices. Both liquid/liquid and gas/liquid element types are 
discussed. Unlike doublet and triplet elements (circular orifices 
only) were investigated for the liquid/liquid case while concentric 
tube elements were investigated for the gas/liquid case. It is 
concluded that noncircular shape can be employed to significant 
advantage in injector design for liquid rocket engines. 
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SYMBOLS 


2 

A = area, in. 

AR = orifice aspect ratio, x^/y^ 

B, b = larger dimension for noncircular orifice, inches 

c* = characteristic velocity, ft/sec 

CCTE = circular concentric tube element 

0^ = specific heat, Btu/lbm-R 

D, d = diameter or characteristic length, inches 

D = mass median diameter, usually microns 

D h = hydraulic diameter, inches 

= mixture ratio uniformity parameter, Eq. F-4 
f( ) = function of ( ) 

f h L = ^ eat loss correction factor, Eq. 3-16 
GCC * gas/liquid circular concentric element code 
GCR = gas/liquid rectangular concentric element code 
G/L - gas/liquid 

GST = gas/liquid multi showerhead triplet element code 

2 

h = heat transfer coefficient, Btu/in. -sec-R 

2 

K = general constant, or bulk compressibility, lbm/ft-sec 

K, k « thermal conductivity, Btu/in. -sec-R 

L = orifice length, or general dimens ion- length 

2 

M = momentum flux, pu , or general dimension-mass 

mff = mass fraction flux 

MR = mixture ratio, w /w- 

o f 
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P, Ap = pressure, psia, and delta pressure, psid 

Pr = Prandtl number 

• 2 
Q = heat flux, Btu/in. -sec 

R = centerpost recess depth, inches 

Re, Rey = Reynolds number 

RCTE = rectangular concentric tube element 

S = wall thickness, inches 

St - Stanton number 

STE = showerhead triplet element 

T = temperature, R, or general dimens ion- time 

U, u = velocity, ft/sec 

U.D. = unlike doublet code 

U. T. = unlike triplet code 

V, v = velocity, ft/sec 

w » smaller dimension for noncircular orifice, inches 


w = flowrate, lbm/sec 

We = Weber number 

x = larger dimension for rectangles, inches 

y = smaller dimension for rectangles, inches 


Greek 

e = contraction ratio 

c 

n = characteristic velocity efficiency 

c* 

U = viscosity, lbm/ft-sec, or microns 
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3 

p = density, lbm/ft 

0 = surface tension, lbf/ft 

2 2 2 

<J> = centerline momentum ratio P^^/P^V d 

Subscripts 


0 

1 

2 

3 

b 


B.C. 


C 


C.F. 


E 

f, F 

g, G 


H.F. 


i 

i 

L 

mix 


o 

opt 


= stagnation properties 

= liquid orifice dimensions, concentric element - smaller orifice, doublets 
outer orifices, triplets 

= centerpost dimensions, concentric element - larger orifice, doublets - 
center orifice, triplets 

= gas port dimensions, concentric element 

= bulk properties 

= coolant bulk 

= coolant 

= cold flow 

= environment properties 
= fuel side properties 
= gas side properties 
= hot fire 

= general index for i quantity 
= jet properties 
= liquid side properties 
= mixing process limited 
= oxidizer side properties 
= optimum condition 
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s = static properties 

sat = saturation properties 

T = total or stagnation properties, also throat conditions 

vap = vaporization process limited 

w = wall properties 

wc = coolant side wall 

wg = gas side wall 
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1.0 SUMMARY 


This report contains the results of a 27-month applied research program. This 
program constituted the final two phases (Phases III and IV) of an overall 4-year 
program. The objectives of Phases III and IV were: (1) to determine the effect of 

shape upon the mixing and atomization characteristics of concentric tube injector 
elements of rectangular shape for gas/liquid propellant system applications, (2) 
to characterize the atomization and mixing performance of liquid/liquid unlike 
triplet elements, and (3) to investigate, parametrically, the effect of area ratio 
and orifice aspect ratio upon the mixing characteristics of liquid/liquid unlike- 
doublet elements with rectangular orifices. (The results of Phases I and II of 
this program may be found in Ref. 1 and 2, respectively, and a comprehensive over- 
all program summary may be found in Ref. 3.) 

The chronology of Phases III and IV effort is presented in Fig. 1-1. Phase III 
was divided into three tasks (all devoted to gas/liquid studies: Task I. Design 

of Gas/Liquid Model Hardware and Definition of the Gas/Liquid Cold-Flow Test Plan; 
Task II. Gas/Liquid Cold-Flow Mixing and Atomization Experimentation; and Task III. 
Design and Fabrication of a Single-element Thrust Chamber for Gas/Liquid Hot-Fire 
Testing Planned for Phase IV. 

Phase IV was divided into five tasks (both gas/liquid and liquid/liquid studies): 
Task I. Hot-Fire Evaluation of the Gas/Liquid Elements Characterized in Phase III 
Cold Flow; Task IA. Evaluation of the Mixing and Atomization Characteristics of 
Unlike Triplets (liquid/liquid cold flow); Task II. Additional Gas/Uquid Cold- 
Flow Mixing Tests; Task III. Evaluation of the Mixing Characteristics of Rectangular 
Unlike Doublets (liquid/liquid cold flow); and Task IV. Final Analysis of results 
and Final Report Preparation. 

1.1 GAS/LIQUID RESULTS 


Concentric tube injector elements of rectangular cross section were investigated 
employing cold-flow techniques. Elements of different aspect ratio, xj/ yi (width to 
height of the central liquid jet), were compared at a fixed set of nominal operating 
conditions. Additional comparisons made included the effects of element size, 
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Figure 1-1. Noncircular Orifice, Phases III and IV Program Schedule 




distribution of gas annuius area, gas-to-liquid area ratio, and center -post recess. 
Results obtained with the rectangular elements were compared with those for a cir- 
cular element of equivalent injection areas. 


Results showed that the mixing uniformity improved markedly with increased aspect 
ratio. For zero center-post recess, the value of E^ increased from 47 percent at 
aspect ratio unity (the circular element) to 65.6 percent at aspect ratio 6 (for 
the nominal conditions) . The improvements provided by aspect ratio were even more 
striking with recessed postsj an of 85 percent was obtained with the aspect 
ratio 6 element at a relative recess, R/y^^ = 1.0.* 

Element size was found to have no effect upon mixing when the flow per unit area 
was maintained. 

Variation of gas port aspect ratio, x^/y^ (g as annulus area distribution about the 
center-post) showed that the characteristics of the flowfield could be altered 
significantly. It was concluded that high levels of mixing could be achieved em- 
ploying this parameter as an independent variable. 

A factor of 2 reduction of the gas-to-liquid area ratio (i.e., an increase in the 
gas-to-liquid velocity ratio by a factor of 2) produced an increase in of 15 
percentage points (60 to 75 percent) . Similar performance trends were noted for 
the circular element. 

It was concluded that the most detrimental parameter to concentric tube mixing is 
high liquid velocity (i.e., high mixture ratio at fixed area ratio). Liquid 
velocities above 10 to 20 ft/sec at mixture ratio 6:1 produce significant reduc- 
tions in E . 

m 


*Post recess also improved performance of the circular element. However, the im- 
provements were greater for the rectangular elements. 
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The effect of liquid jet aspect ratio, x^/y = AR, upon mass median droplet di- 
ameter (frozen wax atomization technique) may be summarized as follows: 


_D 

y i 


AR 2 

3 + AR 



turbulent jets 


where 


K = f (center-post recess) = 0.31 for 2 ero recess 
D = mass median droplet diameter 
y i = height of liquid jet (y^ < x^) 

AR = Xl / yi 


The aspect ratio was varied in such a manner that the liquid port area was main- 
tained constant (i.e., x^/y^ wa.s varied with x^y^ = const). Under this restric- 
tion, the dropsize varied in the following manner: 


K 

D = KV _ 

v£r~ 


AR 


3 + AR 


where A^ - liquid area - const. This relationship shows dropsize increasing with 
aspect ratio to a value of AR~3, and then decreasing with further increase in 
aspect ratio. 

The effect of a reduction of element size, y^ , by 41 percent produced a 70 per- 
cent reduction in dropsize at zero recess. This suggests that the variation of 
dropsize with characteristic element dimension is a complex phenomenon and re- 
quires more than the relative droplet diameter, D/y^, for description. 

The redistribution of gas annulus area from uniform gap to nonuniform gap around 
the central liquid post produced a reduction in dropsize. It is expected that 
this is not a general trend but depends upon the degree and the direction of the 
redistribution of annular gap. 
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Gas/liquid hot fire results (4-inch long, e c = 2.3 single-element chamber) were 
shown to be nearly mixing limited with the vaporization efficiency being approx- 
imately unity under the design and operating conditions tested. The extremely 
high vaporization efficiency are presently attributed to secondary droplet breakup. 
This effect was not included in predicted vaporization efficiencies. The hot 
values are, therefore, higher than those calculated. The trends of the hot-fire 
results confirmed most of the mixing trends suggested by the cold-flow results. 

It was generally concluded that the rectangular concentric tube element could 
most advantageously be employed to improve mixing efficiency for high mixture 
ratio, and/or low gas momentum injector applications. The rectangular elements 
showed no particular atomization advantage. This result suggests that the atomiza- 
tion process is strictly governed by the energetics of the shearing process and is 
little effected by the relative shape of the flow field. 

1.2 LIQUID/LIQUID STUDIES 


Triplet mixing results (circular orifices) showed that optimized for a given 
element at a value of centerline momentum ratio of unity for diameter ratio 
dj /^2 = 1 -°* Cdj_ is the outer orifice diameter and d 2 is the central orifice 
diameter.) Optimum mixing occurred at 0 = 0.7 for d 1 /d 2 < 1.0. The parameter 0 
has the following definition for triplets: 


p i V d i 

0 = L 1 1 

P 2 V 2 2 d/2 


centerline momentum ratio 


It was also noted that the value of E m dropped as the value of dj/d.^ was de- 
creased from unity. °P t 

Only laminar jet atomization data were correlated for the triplet elements. The 
data were shown to follow a functional relationship of the form: 


D_ 

D. 

J 


f (Weber number x — 


) 
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where 


D 


mass median droplet diameter 
given orifice diameter 


p V 2 Dj 

Weber number = _ 

a 

L/D = orifice length-to-diameter ratio 

The functional relationship between dropsize and physical parameters was found 
to be: 


D = 


K V .V-. 77 / 


(L/D) 


-.39 


.62 | 


where K is a constant. It was concluded that the relatively strong influence of 
L/D was due to its affect on jet velocity profile development as suggested by the 
results of Zajac (Ref. 9). 


Results of the investigation of the mixing characteristics of unlike-doublet 
elements with rectangular orifices suggested the following design criteria for 
optimum mixing: 


and 


0 = 1 




8 


where 


0 


P f V f 


P v 
o o 


b 


b 


f 


o 
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with the following definition of geometrical parameters for a rectangular unlike 
doublet: 




1 

W 

T 


These results are based upon cold flow and are not directly indicative of the 
impingement characteristics of hypergolic reactants that can be dominated by 
reactive stream separation or "blowapart 
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2.0 INTRODUCTION 


The study of the influence of orifice and element shape upon rocket engine injec- 
tor performance was initiated under the Nonciruclar Orifice Program in June of 
1969 under NASA funding from the Manned Spacecraft Center (now the Lyndon B. 
Johnson Spacecraft Center), Houston, Texas. In Phase I of this program (Ref. 1) 
the flow characteristics of individual orifices of noncircular shape and the mix- 
ing and atomization characteristics of unlike-doublet injector elements with non- 
circular orifices were investigated. The scope of Phase I was limited to liquid/ 
liquid propellant systems, primarily of the storable, hypergolic variety. It was 
discovered that the noncircular shape offered greater design flexibility for the 
unlike-doublet element by providing high levels of mixing uniformity over a broad 
range of propellant density and mixture ratio. This flexibility had, heretofore, 
been unavailable with circular orifices. 

The scope of the program was expanded in Phase II to include gas/liquid propel- 
lant combinations (Ref. 2) with emphasis upon the L0X/GH 2 system. Phase II con- 
sisted of an analytical screening analysis of candidate gas/liquid elements to 
select interesting elements for further research and a preliminary experimental 
investigation of a concentric tube injector element of rectangular shape. The 
screening study indicated that the concentric tube element type had the greatest 
potential for improvement with the adaptation of noncircular shape, and had the 
largest data base for comparison. Preliminary cold-flow experimental results sug- 
gested that improvements in both atomization and mixing could be achieved with 
rectangular shape. The atomization results were quite straightforward. However, 
the mixing uniformity trends with physical parameters were unclear. 

Results obtained under both Phases I and II stimulated interest in further in- 
vestigation of noncircular orifices for both liquid/liquid and gas/liquid propel- 
lant systems. Rocketdyne was subsequently funded to conduct a Phase III and 
Phase IV effort. The scope of the additional research was formulated to include 
extensive cold-flow mixing and atomization studies for gas/liquid concentric tube 
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elements, hot-fire investigation of gas/ liquid concentric tube elements, and 
further cold-flow mixing and atomization studies with (liquid/ liquid) unlike 
doublets and triplets. 


This volume is a presentation of the results of the Phase III and IV portion of 
the Noncircular Orifice Program. A comprehensive summary of the entire program 
is contained in Ref. 3. 
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3.0 GAS/LIQUID STUDIES 


The gas/liquid studies reported herein are composed primarily of cold-flow mixing, 
cold-flow atomization, and hot-fire experimental investigations of the effect of 
geometrical parameters upon the performance of rectangular concentric tube injec- 
tor elements. (A small amount of cold-flow data was obtained for another element 
type, a showerhead triplet. Reporting of the results of that element is limited 
to a discussion at the end of the gas/liquid section to avoid confusion of the 
concentric tube presentation with another type of geometry.) To avoid unnecessary 
repetition of certain lengthy terms, a list of abbreviations to be employed 
throughout this report is presented in Table 3-1. 


TABLE 3-1. LIST OF ABBREVIATIONS FOR GAS/LIQUID 
INJECTOR ELEMENT STUDIES 


Term 

Abbreviation 

Gas/Liquid 

G/L 

Rectangular Concentric Tube Element 

RCTE 

Circular Concentric Tube Element 

CCTE 

Showerhead Triplet Element 

STE 


3.1 TECHNICAL APPROACH 

The technical approach adopted for the G/L study is outlined in block diagram 
form in Fig. 3-1. The objective of the program was to determine experimentally 
the influence of element geometry upon injector performance employing rectangular 
concentric tube elements (RCTE), and, as a result, to answer the question, 

" Exactly what can a noncircular element accomplish that cannot he accomplished 
by a circular element?" To attack this problem, the concept of "performance" 
had to be broken down into fundamental components; one related to the quality 
of mixing provided by a given element, and one related to the degree to which 
that element is capable of atomizing the liquid propellant. This was accomplished 
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Figure 3-1. Technical Approach for the Characterization of Gas/ 
Liquid Rectangular Concentric Tube Injector Elements 
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by assuming that thrust chamber c* efficiency is approximated, to the first order, 
by the product of a mixing-limited and a vaporization-limited c* efficiency 
(Fig. 3-1 and Appendix F). Under this assumption, the performance problem is 
greatly simplified and can be investigated along two essentially independent paths. 

The next step in the approach was to single out those geometrical and operational 
parameters that, (1) most influence the mixing and atomization processes, and 
(2) could be varied over significant ranges employing practical experimental tech- 
niques within the scope of the program. This was accomplished by the application 
of dimensional analysis. 

With the selected parameters as guidelines, a set of model hardware and an exper- 
imental approach were formulated. The model hardware was designed to allow vari- 
ation of the selected geometrical parameters. The experimental plan included 
independent cold-flow mixing, cold-flow atomization, and hot-fire experimentation. 

Primary element comparisons and design criteria are obtained from the cold-flow 
results. The hot-fire experimentation was performed to ensure that the trends 
indicated by the cold-flow results were essentially correct. 

To be able to compare cold-flow and hot-fire performance trends, the cold-flow 
results were converted to equivalent mixing-limited and vaporization-limited c* 
efficiencies (see Appendix F). The product of these efficiencies was termed 
c anc * was com P are d to the hot-fire efficiency, q „ . 

It was not the objective of this study to "predict" performance. However, the 
hot-fire trends must serve as a yardstick for the credibility of the cold-flow 
design criteria. 

A CCTE was carried along throughout the program to serve as a basis of comparison 
for the RCTE results. It should be emphasized that the comparisons between cir- 
cular and rectangular elements made in this report are not intended to imply that 
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circular elements are either "worse” or "better" than noncircular over all ranges 
of injector design. These results should be viewed in the light of changes in 
performance (at the conditions and over the ranges stipulated in this document) 
produced by variations of shape, alone. 

3.2 DEFINITION OF SIGNIFICANT PARAMETERS AND 
TEST PLAN FORMULATION 

3.2.1 Dimensional Analysis 


As could be envisioned, intuitively, the mixing and atomization processes associ- 
ated with concentric tube injector elements are highly complex. To single out 
the individual parameters that should most influence these processes, the method 
of dimensional analysis was applied to the problem. This method, attributed to 
Buckingham, stipulates that if there are M physical parameters involving N dimen- 
sions, then there are M-N dimensionless groups that will completely characterize 
the physical problem. 

The first task in the development of a dimensional analysis is to identify all of 
the parameters that are involved in the problem. Great care must be taken to en- 
sure that all of the parameters are included. At this point, a discussion of the 
flowfield and the physics associated with concentric- tube element mixing and 
atomization is appropriate. 

Injector face geometries for both a CCTE and a RCTE are shown in Fig. 3-2, and a 
cross section of a concentric tube element with physical processes identified is 
presented in Fig. 3-3. The standard concentric -tube element (often called coax- 
ial element) is composed of a single tube within a larger orifice. Liquid, usu- 
ally oxidizer, is injected through this central tube while gas, usually fuel, is 
injected through the annulus between the central post and the orifice (Fig. 3-3). 
Atomization of the liquid and subsequent mixing of the gas and liquid is produced 
by shear between the gas and liquid. The bulk of the energy required for the mix- 
ing and atomization processes is provided in the gas component in the form of 
velocity head. 


R-9270 

14 



(a) CIRCULAR 



(b) RECTANGULAR 


Figure 3-2. Typical Face Geometry for Circular 
and Rectangular Concentric Tube 
Elements 


R-9270 
15 



R-9270 


ENVIRONMENTAL GAS 


MIXING BETWEEN PRIMARY AND 
ENVIRONMENTAL GAS 



Figure 3-3. 


Typical Concentric Tube Element Flowfield 



In many cases, the performance of the elements may be improved by recessing the 
central tube (Fig. 3-3) to allow more efficient momentum (or energy) exchange 
between the two fluids. There are situations, however, in which recess can be 
detrimental either from a performance or a hardware compatibility standpoint. 

With the physical processes in mind, it is very important to note that there is 
a significant difference between a cold-flow model of these processes and the 
actual operation of an element with chemically reactive propellants. In the cold- 
flow situation, there is no vaporization of the liquid component and there is no 
energy release due to combustion. Thus, cold-flow results should not be expected 
to "predict" levels of performance exactly. However, cold-flow trends in perform- 
ance with the geometrical and operational variables should be representative of 
actual trends. These modeling restrictions are not as important in the cold-flow 
modeling of liquid/liquid injector processes due to the fact that the majority of 
the mixing and atomization of liquid/liquid propellants occurs before significant 
combustion interaction. 

In addition to the primary mixing of the injected liquid and gas components, there 
is a competitive mixing of the environmental gas with the injected gas. This sec- 
ondary mixing is detrimental to the performance of the element as it severely de- 
grades the momentum of the primary gas jet. The importance of the secondary mixing 
is diminished, however, with increased central tube recess. 

The variables that were selected as representative of the concentric-tube element 
performance are listed in Table 3-2, along with their basic physical dimensions. 
(Note must be taken that many variables have been eliminated from the analysis. 

For example, the compressibility of the liquid has been assumed to be constant, 
gravitational forces have been neglected, and all temperatures, specific heats 
and, thus heat transfer, have been neglected.) There are 19 variables considered, 
having a total of three dimensions. According to the Buckingham II theorem (see, 
e.g. , Ref. 4), with M = 19 variables and N = 3 dimensions, there are 19 - 3 = 16 
dimensionless groups that will completely characterize the problem. 
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TABLE 3-2. PHYSICAL VARIABLES AND THEIR DIMENSIONS 


t 

Variable 

Liquid 

Gas 

Environment 

Dimensions* 

! 

Geometry 

Xy X 2 , X 3 , y 1# y 2 , Vy R 

L 

Density 

P L 

P G 

P E 

ML' 3 

Velocity 

V L 

V G 

V E 

LT" 1 

Compressibility 


k g 

k e 

ML" 1 T' 2 

| — ■ 

Viscosity 

y L 

y G 

V E 

ml" 1 t -1 

Surface Tension 

°L 



mt" 2 


M = mass 
L = length 
T = time 
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The results of the dimensional analysis are presented in Table 3-3. These groups 
are the fundamental groups. More meaningful combinations can be obtained by com- 
bining two or more groups and replacing one of the variables employed in the new 
group. For example, the mixing of the primary gas with the environment is a free 
turbulent mixing process and is little dependent upon the viscosity of the gases 
(Reynolds No.). A better group may be (p Q V G - P £ V £ )/p G V G to replace V £ /V G . 
However, no attempt has been made in this report to regroup the variables, a prion. 
Rather, the quantities which appear in the groups that could be varied were studied, 
and the results used to suggest which dimensionless groups would be most appropriate. 

TABLE 3-3. DIMENSIONLESS GROUPS REQUIRED FOR THE DESCRIPTION OF GAS/ 

LIQUID CONCENTRIC-TUBE ELEMENT MIXING AND ATOMIZATION 


Group 

No. 

Group 

Name and Type 

1 

x l^l ( As P ect R atio ) 


2 

x 2 /y 2 


3 

V y 3 

Geometry 

4 



5 

y 3 /y 2 


6 

R/yj (Relative Recess) 


7 

P CT / Pr 

Density Ratio 


g L 

8 

p e/ p g 

9 

V V T 



G L 

Velocity Ratio 

10 

V E^ V G 

11 

P G V G ^ y 3~ y 2^ P G 


12 

P L V L (y l )/u L 

Reynolds Number 

13 

P E f VV y 3 /U E 


14 

V' k g /p g 

Mach Number 

15 

! 

V E^ K E^ P E 

mm 

p G <W 2 y l /0 L 

Weber Number 
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It is obvious upon inspection of the results shown in Table 3-3 that the task of 
investigating the concentric-tube flowfield is formidable. Therefore, the next 
task in the formulation of the experimental approach was to significantly reduce 
the number of variables for further study. 

Since most of the basic results to be generated under this contract would come 
from cold flow, it was natural, first, to examine which of the parameters involved 
in the problem could be varied over a significant range and yet maintain the study 
within the scope of the contract. 

Realistically, the only parameters that could be varied included the geometry, the 
gas density, the gas velocity, and the liquid velocity. Therefore, these param- 
eters were selected as the primary experimental variables. Since all the dimen- 
sionless groups could not be varied independently, it was decided that presentation 
of the results in terms of any of these groups would be inappropriate unless sound 
justification could be presented. Thus, only in the case of certain of the geo- 
metric groups are results presented in dimensionless form. However, where possible, 
the results are employed to suggest which terms should be investigated in more 
detail . 

In addition to these restrictions, it was also decided that the scope of the pro- 
gram did not allow an independent variation of all of the geometrical groups. The 
rationale employed to limit the experimental scope and the basis of the test plan 
are presented in the next section. 

3.2.2 Test Plan Formulation 


In formulating the test plan, the most important consideration was the original 
objective of the study: determine the influence of shape upon the performance of 

concentric-tube injector elements. Therefore, the question that one should ask 
is: "Exactly what does a rectangular concentric-tube element accomplish that a 

circular element cannot?" 
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Study of Table 3-3 will show that all of the dimensionless groups can be varied 
equally well by either a CCTE or a RCTE except groups No. 1, 2, and 3. It is clear, 
upon reference to Fig. 3-2 that these groups are available only to RCTE's. Thus, 
these three groups embody the only possible advantages or disadvantages available 
with noncircular (rectangular) shape. 

The experimental plan, therefore, was limited to a parametric variation of these 
geometrical groups at a nominal baseline set of operating conditions. Some vari- 
ations about the nominal conditions were selected to show relative sensitivity of 
the various elements. However, the restriction to baseline conditions was neces- 
sary, again, to meet scope. 

Operating conditions selected as baseline are those of an 800-psia chamber pressure, 
liquid oxygen, gaseous hydrogen (LOX/GH^) rocket engine with contraction ratio ap- 
proximately 2.3:1, and nominal flow per element of 0.43 Ib/sec. A list of the base- 
line conditions is presented in Table 3-4. 

TABLE 3-4. BASELINE OPERATIONAL PARAMETERS 


Parameter 

Symbol 

Nominal Value 

Chamber Pressure 

P 

c 

800 psia 

Propel lant Properties 


lox/gh 2 

Mixture Ratio 

MR 

6 

Gas Density 

P G 

0.27 lbm/ft 3 

Liquid Density 

p L 

70 lbm/ft 3 

Chamber Contraction Ratio 

e 

c 

~2. 3 

Flow per Element 

Wp 

0.432 lbm/sec 

Fuel Injection Velocity (gas) 

V G 

~1000 ft/sec 

Oxidizer Injection Velocity (liquid) 

V L 

^50 ft/sec 
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The set of concentric-tube elements selected for experimental evaluation is shown 
in Fig. 3-4 to relative scale. (Letters and numbers appearing under each element 
are element code numbers.) 

The element comparison was accomplished in two stages: (1) comparison between a 

CCTE and "baseline" RCTE's, and (2) a comparison of one baseline RCTE with three 
variations of geometrical parameters about this baseline (Fig. 3-4). 

The primary geometrical variable selected for comparison is the liquid orifice 
dimension ratio, Xj/yj termed the element aspect ratio (AR) , (Fig. 3-2). Values 
of 1.5, 3.0, and 6.0 were chosen for this ratio. Certain ground rules were estab- 
lished for the primary baseline comparison; all baseline liquid areas were equal, 
all baseline gas annulus areas were equal, gas annulus gap was uniform, and all 
central tube wall thicknesses were equal. This was stipulated to ensure that all 
baseline elements would produce respectively equal gas velocities and liquid veloc- 
ities at the baseline operating conditions. The equivalence of post tube-wall 
thickness was a fabrication consideration. 

Under these ground rules, all the geometrical dimensionless groups become directly 
related to the aspect ratio, X^/y^ 

The central post aspect ratio of 3.0 was selected for further parametric study. 
These additional elements are labeled "variations about AR = 3," in Fig. 3-4. One 
of these elements, GCR-7 is a direct scale model of the baseline AR = 3 element, 
reduced by a factor of 1//2" . This reduces all of the areas by a factor of 
two. The other two variations (GCR-4 and GCR-3) employ the exact center post of 
GCR-2 (basic AR = 3 element) and have, respectively, one-half the gas annulus area 
(uniform gap) and a nonuniform gas gap (same gas area as GCR-2, Fig. 3-4). 

These latter three variations were designed to test the effects of element size, 
velocity ratio, and nonuniform gas annulus distribution. 
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GCR-7 QCR-4 GCR-3 

SCALED DOWN 1/2 GAS AREA NONUN I FORM GAP 

BY 1//2 


Figure 3~4. Rectangular Concentric Tube Element Comparisons 
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A complete list of all of the elements, their dimensions, and the objective of 
each is presented as Table 3-5. Two additional elements (GST-1 and GST-2) appear 
in Table 3-S. These are the "showerhead triplet elements" and are discussed in a 
separate section at the end of major section 3.0. The code number meaning is ex- 
plained in Table 3-5. 

In the next section, the hardware designed to implement the test plan is described 
followed by the presentation of the experimental results. The test plan described 
in this section has been general in nature. The exact test plan and method of 
presentation are detailed in Section 3.4. 

3.3 TEST HARDWARE 


The discussion of the test hardware is divided into two sections: injectors and 

^°t~fire thrust chamber. The injectors were designed for use in both the cold- 
flow and the hot-fire portions of the program; therefore, there is further dis- 
cussion of injectors directly in the thrust chamber section. 

3. 3. 1 Injectors 

Detailed injector face geometry descriptions and rationale for geometry selection 
have been presented and are summarized in Fig. 3-4 and Table 3-5. Only the hard- 
ware aspect of the injectors is presented in this section. 

As mentioned, the injector models were designed to serve both the cold-flow and 
the hot-fire tasks. Each injector is composed of four pieces (Fig. 3-5)--a body, 
a clamp (or retaining ring), a post, and a port. The body and the clamp were 
employed for all of the models, only the port and/or post were changed to achieve 
a different configuration. A photograph of the individual components for a typ- 
ical element is presented in Fig. 3-6 (element GCR-5 is shown; Table 3-5). 

The body holds the entire assembly together and serves as a manifold for the gas- 
eous component. The post is held into the body by means of the retaining clamp 
ring that bolts to the body. The post is maintained "concentric" with the port 


R-9270 

24 



R-9270 


TABLE 3-5. GAS/LIQUID ELEMENT GEOMETRY 
(See Fig. 3-2 and 3-4) 


Element * 
Code 
No. 

X 1 

(in.) 

X 2 

(in.) 

X 3 

(in.) 

y i 

(in.) 

y 2 

(in.) 

y 3 

(in.) 

Aspect 

Ratio 

V y i 

Gas 

Area 

(x 3 / y 3 - 

x 2 /y 2 ) 

Liquid 

Area 

x i /y i 

Element Description 

GCC-l 

0.160 

0.200 

0.305 

0.160 

0.200 

0.305 

1,0 

0.0416 

0.0201 

Circular concentric tube from 
Phase II 

GCC- 2 

0.1604 

0 . 2005 

0.292 

0.1604 

0,2005 

0.292 

1.0 

0.0354 

0.0202 

Circular concentric tube from 
Phase III 

GCR- 1 

0.346 

0.386 

0.450 

0.058 

0.098 

0.162 

5.97 

0.0351 

0.0201 

Rectangular concentric tube from 
Phase II, AR=6 

GCR-2 

0.245 

0.286 

0.3592 

0,0825 

0.1211 

0.1936 

2,97 

0.0349 

0.0202 

Baseline RCTE, uniform gap, AR=3 

CCR-3 

0.245 

0.286 

0.4039 

0.0825 

0.1211 

0.1721 

2.97 

0.0349 

0.0202 

Nonuniform gas gap, AR=3 

GCR-4 

0.245 

0.286 

0.3249 

0.0825 

0,1211 

0.1615 

2.97 

0.0178 

0.0202 

1/2 baseline gas area, AR=3 

GCR-S 

0.350 

0.3876 

0.4518 

0.0595 

0.0981 

0.1625 

5.88 

0.0354 

0.0208 

Baseline RCTE, uniform gap, AR=6 

CCR-6 

0.173S 

0.2122 

0.2922 : 

0.1165 

0.1611 

0.2333 

1.49 

0.0340 

0,0202 

Baseline RCTE, uniform gap, AR=1.5 

GCR- 7 

0.1750 

0.2020 

0.2535 

0.058 

0.0851 

0.1379 

3.02 

0.0178 

0..0101 

1//2 scale model of GCR-2 

GCR- 8 

0.2450 

0.286 

0.3590 

0,0825 

0.1211 

0.1950 

2.97 

0.0354 

0.0202 

Remake of GCR-2 with copper gas 
port for hot fire 

GCR- 9 

0.1735 

0.2122 

0.2903 

0.1165 

0.1611 

0.2350 

1.49 

0.0340 

0.0202 

Remake of GCR-6 with copper gas 
port for hot fire 

GCR-10 

0.350 

X J 

(in.) 

0.3876 

0.450 

0.0595 

i y i 

(in.) 

0,0981 

D. 

J 

(in.) 

0.1620 

5,88 

0.0349 

0.0208 

Remake of GCR* 5 with copper gas 
port for hot fire 

GST-1 

0.320 

--- 

— 

0.055 

0.020 


5,82 

0.0352 

0,0088 

Showerhead triplet with 28 each 
0.020 in. dia. orifices 

GST- 2 

0.320 



0.055 

0.030 


5.82 

0.0352 

0.0088 

Showerhead triplet with 12 each 
0.030 in, dia. orifices 


# Code No* Definition 

GCC = Gas/Liquid Concentric Tube - Circular 
GCR = Gas/Liquid Concentric Tube - Rectangular 
GST = Gas/Liquid Showerhead Triplet 
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1DZ36-6/21/71-S1C 


Figure 3-6. Injector Components (Element GCR-5 Shown) 



means of the alignment tabs, which may be seen in Fig. 3-5 and 3-6 toward the 
downstream end of the post. These tabs are relatively far from the gas annulus 
in a low-velocity region and did not pose a significant wake problem in the gas 
flow. The port is bolted to the body. 

Sealing of the liquid, which is injected straight through the post, is accomplished 
by means of standard AN type B-nut fitting (not shown). Sealing of the gas is 
accomplished by two O-ring seals (Fig. 3-5). The gas is injected through two tubes 
that were brazed to the gas inlet orifices in the body. (Only one tube is shown in 
Fig. 3-5 due to the manner in which the section was taken.) 

Center-post recess capability was provided by placing spacers between the body and 
the ring on the posts. 

The liquid orifices were provided with rounded inlets to avoid flow separation and 
cavitation and were designed with a length to minimum dimension (yj) ratio of 10:1. 
The gas annulus gaps were designed with a length-to-gap height ratio of 5:1 at zero 
recess conditions. 

The i ectangular orifices in the ports and posts were produced by the electrical 
discharge machining (EDM) method. This process yields an orifice with very accur- 
ate dimensions but with a rather rough (sand-blasted- like) surface. 

All components of the injectors were made of series 321 stainless steel except the 
port pieces employed in the hot-fire experiments. These were made from OFHC copper 
for improved heat transfer capability. 

Front and back views of the injector assembly are shown in Fig. 3-7 and 3-8, re- 
spectively (element GCR-2 shown). The two gas inlet orifices can be seen in 
Fig. 3-8 before brazing of the inlet tubes. 

A comparative photograph of the three baseline RCTE’s (Table 3-5) is presented 
in Fig. 3-9, wherein only the ports and posts are shown. 
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Figure 3-7. Injector Assembly, Front View (Element GCR-2 Shown) 
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Figure 3-8. Injector Assembly, Rear View (Element GCR-2 Shown) 
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Figure 3-9. Comparison of the Three Baseline Rectangular Concentric-Tube Elements 


A face detail of element GCR-2 is shown in Fig. 3-10, at high magnification, to 
illustrate the dimensional control which was achieved with the EDM process and 
the concentricity between the port and post provided by the alignment tabs on 
the post. 

The showerhead triplet injector models (GST-1 and GST-2) employ the same body and 
retaining clamp as the RCTE's and the CCTE. However, instead of a liquid post, a 
manifolding piece is used to introduce the liquid to the injection orifices. A 
cross section of this element is shown in Fig. 3-11 and a photograph of the face 
of GST-1 is shown in Fig. 3-12. The liquid is injected through many small shower- 
head orifices (circular), while the gas is injected through two orifices directed 
at the liquid jets with an included angle between the gas jet orifice centerlines 
of 60 degrees. The gas orifices are rectangular with aspect ratio of 5.82:1. 

Gas is supplied to these orifices through the same manifold employed with concen- 
tric tube elements. There is an additional 0-ring seal between the liquid and 
the gas at the tip of the center-post manifold (Fig. 3-11). 

All injector models were designed to adapt to the pressurized cold-flow mixing 
facility (Appendix A), the cold-flow atomization facility (Appendix B), and the 
hot-fire thrust chamber (Section 3.3.2). 

For documentation, all Phase III gas/liquid element hardware is keyed to Rocket- 
dyne drawing numbers in Table 3-6. It will be noted that elements GCC-1 and 
GCR-1 do not appear. These elements were employed in Phase II and were of dif- 
ferent design (Ref. 2). 

3.3.2 Hot-Fire Thrust Chamber 


The discussion of the single- element , hot-fire thrust chamber design is presented 
in two sections--a summary of the heat transfer analysis for the chamber followed 
by a discussion of the chamber design itself. 
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5AA32- 8/20/71 -SI 


Figure 3-10. Element Face Detail (Element GCR-2 Shown) 



Figure 3-11. Single-Element Injector Assembly (Element GST-2 Shown) 
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5AA34- 8/2/71- SI 


Figure 3-12. Multishowerhead Triplet Face Pattern (Element GST-1 Shown) 



TABLE 3-6. PHASE III ELEMENT HARDWARE DRAWING NUMBER KEY 



Drawing Numbers 

■C JL viiitll L 

Code 

Body 

Clamp 

Port 

Post 

GCC-2 

Ric; 

2928 

RIB 

2930 

R1B2960 

R1C2961 

GCR-2 





R1B2927 

R1C2929 

GCR-3 





R1B29S4 

R1C2929 

GCR-4 





R1B2955 

R1C2929 

GCR-5 





R1B2935 

R1C2933 

GCR-6 





R1B2934 

R1C2931 

GCR-7 





R1B2957 

R1C2958 

GCR-8 





R1B2927-3 

R1C2929 

GCR-9 





R1B2934-3 

R1C2933 

GCR-10 





R1B2935-3 

R1C2933 

GST-1 





R1B2937 

R1C2938 

GST- 2 


r 



R1B2948 

R1C2938 
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3. 3. 2.1 Chamber Heat Transfer Analysis . The chamber heat transfer analysis was 
based upon the nominal operating conditions listed in Table 3-4. With the assump- 
tion of 100-percent c* efficiency, these conditions yield the chamber design param- 
eters listed in Table 3-7. 

A chamber length (injector to throat) of 4 inches was selected based upon a pre- 
liminary droplet vaporization calculation that indicated that vaporization would 
be incomplete at the throat. Thus, performance variations due to dropsize varia- 
tions were expected to be clearly evident within the accuracy of the experimental 
apparatus . 

A heat transfer analysis of the chamber was conducted with Rocketdyne's integral- 
momentum turbulent boundary layer computer program. The chamber wall temperature 
was assumed to be a constant 1560 R. The results of that analysis are presented 
in Table 3-8. 

A copper chamber liner of the type shown in Fig. 3-13 with axial coolant channels 
was selected for the chamber design. Subcooled water was chosen for the coolant. 

The design of the copper chamber liner was based upon four heat transfer correla- 
tions and/or fundamental equations plus a stress analysis of the liner wall at 
the thinnest section. These design equations are listed in Table 3-9. 

A great deal of judgment and experience was applied to the design in addition to 
the equations listed in Table 3-9. Many of the variables appearing in the equa- 
tions are either known or can be estimated to the first order. 

The following is a brief discussion of the logic employed to reduce the number of 
variables in the design equations. It is arranged in list form for clarity of 
presentation. 
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TABLE 3-7. PRELIMINARY CHAMBER DESIGN PARAMETERS 


Parameter Value 


Maximum Adiabatic Temperature 

6374 R 

Throat Diameter 

0.395 inch 

Chamber Diameter 

0.600 inch 

Pu at Throat 

3.525 lbm/in. 2 -sec 

pu in Chamber 

1.528 lbm/in. 2 -sec 

Throat Reynolds Number 

2.62 x 10 6 

Chamber Reynolds Number 

1.72 x 10 6 

Specific Heat Ratio 


Frozen 

1.202 

Shifting 

1.138 

Specific Heat (cp) 


Frozen 

0.890 Btu/lbm-R 

Shifting 

2.417 Btu/lbm-R 

Prandtl Number 

0.826 


TABLE 3-8. RESULTS OF THE TURBULENT BOUNDARY LAYER ANALYSIS 
(P = 800 psi; T 0 = 6374 R; T w = 1560 R) 


Parameter 

Value 

Throat Heat Flux 

2 

33 to 35 Btu/in. -sec 

Chamber Heat Flux 

2 

19 to 20 Btu/in. -sec 

Gas-Side Film Coefficient 


Throat 

0,00685 Btu/in. 2 -sec-R 

Chamber 

0.00414 Btu/in. -sec-R 
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CHAMBER SECTION 



Figure 3-13. Cross-Sectional View of Hot-Fire Chamber Liner 
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TABLE 3-9. HEAT TRANSFER RELATIONSHIPS EMPLOYED FOR THE DESIGN OF THE COPPER CHAMBER LINER 

(ONE-DIMENSIONAL CONDUCTION ASSUMED) 


Relationship 

Equation* 

Eq. No. 

Reference 

Heat Flow Continuity 

4- ‘v - v - £ 

3-1 

Basic 

Definition 


h c (T wc ‘ T BC } = i 

3-2 


Bum Out Heat Flux 

• 01 ^c‘ T S AT- T BC) . (?) 

\ a /b.o. 

3-3 

Ref. 5 

Coolant Side Heat 
Transfer Coefficient 

* -,05 -.6 

S t = .005 (R e ) (P r ) 

3-4 

Ref. 6 


* k = Conductivity of Liner Wall, BTU/in-sec-°R 
s = Liner Wall Thickness, inch 
* h = Coolant Side Coefficient, BTU/in -sec-°R 
Q/A - Actual Estimated Heat Flux, BTU-in -sec 
(Q/A)B.O. - Burnout Heat FluXj BTU/in -sec. 

V c « Coolant Velocity, Ft/Sec 


T = Gas Side Wall Temp., °R 

r® = Coolant Side Wall Temp., °R 

■T.WC 

a BC = Coolant Bulk Temp., °R 
T qAT = Coolant Saturation Temp., °R 
S a h / Pvc (Stanton) 

R = p vd/y p (Reynolds) 

P » cM/k (Prandtl) 

r p 




1 . 


The value of Q/A maximum was given by the gas-side heat transfer solution 

* 2 
Q/A = 35 Btu/in. -sec. 
x max 

2. The burnout heat flux was computed by multiplying the maximum heat flux 

by an assumed cooling margin. A value of 1.25 was selected for a cooling 

2 

margin; therefore (Q/A)g ^ = 1.25 x 35 = 44 Btu/in. -sec. 

3. With copper liners, it is reasonable to limit the gas-side wall tempera- 
ture to 800 F or 1260 R; therefore T w g = 1260 R. 

4. The saturation temperature for water is shown plotted as a function of 
static pressure in Fig. 3-14. It is desirable to maintain a high satura- 
tion temperature while limiting the pressure to reasonable values. For 
pressures of the order of 600 psia, the saturation temperature is of the 
order of 500 F; therefore Tg^ = 500 F = 960 R. 

5. For added yield strength, at elevated temperature and pressure difference, 
NARloy was selected for the chamber liner material rather than OFHC copper. 
Thus, the thermal conductivity was established, K = 0.005 Btu/in. -sec-R. 

6. The facility for hot-fire experimentation could be expected to supply 
coolant water at temperatures around the ambient temperature. Actually, 
a subsequent analysis, which is presented later in this section, indi- 
cated that a coolant water temperature of the order of 120 F was superior 
to ambient temperatures; thus, = 580 R. 

7. The scale of the liner dictates that the hydraulic diameter of the chan- 
nels should be of the order of 0.1 inch. Reasonable coolant velocities 
could be expected to be roughly 150 ft/sec. These approximate values 
yield a coolant Reynolds number of the order of 2 x 10 . 

8. Around 100 F, the Prandtl number for water is about 4.0;- Pr = 4.0. 
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STATIC PRESSURE, PS I A 

Figure 3-14. Saturation Temperature of Water as a Function of Pressure 
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These approximate values are listed in Table 3-10. Substitution of the values 
from Table 3-10 into the equations in Table 3-9 yields the following set of sim- 
plified equations. The coefficient in Eq. 3-4 was reduced by 60 percent for 
safety margin, and for water was assumed to be 1 Btu/lbm-R. (S is the wall 
thickness in inches.) 


T wc = 1260 - 7000 S 

(3-1) 

h c (T wc - 580) = 35 

(3-2) 

V = 134 

c 

(3-3) 

h = 0.000307 V 

c c 

(3-4) 


Solution of these equations gave the first-order liner design parameters: 

V c = 134 ft/sec 

h c = 0.0411 Btu/in. ^-sec-R 

T wc - 1432 R 

S = 0.0246 inch 

These numbers are, obviously, not physical. The problem rests with allowing the 
wall thickness to vary, arbitrarily, with the rough-order-of-magnitude estimates 
for the other variables. 

If the wall thickness, S, is selected from structural considerations, it can be 
seen that the set of design equations is overdetermined. That is, there are two 
criteria for selection of the coolant velocity; one from burnout considerations, 
and one from coolant-side heat transfer conditions. If the velocity required to 
achieve the coolant-side coefficient is greater than the burnout velocity, then 
the design parameters are acceptable to the first order. 
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TABLE 3-10. APPROXIMATE VALUES EMPLOYED FOR 


"FIRST-TRIAL" CHAMBER LINER DESIGN 


Variable 

Approximate Value 

* 

Q/A 

35 Btu/in. ^-sec 

«W B .o. 

2 

44 Btu/in. -sec 

T 

1260 R 

wg 


T 

960 R 

sat 


K 

0.005 Btu/in. -sec-R 

T 

BC 

580 R 

r 

Re 

2 x 10 s 

Pr 

4 


A structural analysis of the liner wall was conducted and it was concluded that 
a minimum liner wall thickness equal to 0.046 inch was acceptable. Setting S 
equal to 0.046 inch in Eq. 3-1 and solving the set of design equations gave the 
following results: 

T wc = 938 R = 478 F 

h^ = 0.0978 Btu/in. ^-sec-R 

(V ) = 134 burnout 

(V ) = 318 coolant-side coefficient 

Based upon these preliminary results, it was concluded that coolant-side coeffic- 
ients of the order of 0.09 Btu/in. -sec-R and coolant velocities of the order of 
at least 200 ft/sec at the throat must be employed. A chamber liner design was 
configured to meet these criteria. The liner is shown in cross section in 
Fig. 3-13. Eighteen evenly spaced axial coolant channels were selected, having 
a total flow area of 0.0591 in. in the chamber section and 0.04536 in. at the 
throat. A coolant velocity of 200 ft/sec was selected, which gave a total water 
flowrate of approximately 4 lbm/sec and a coolant velocity of 153 ft/sec in the 
chamber section. Wall thickness at the channel roots, S, was set at 0.046 inch. 
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The liner is shown in cross section in Fig. 3-15, along with the heat fluxes, 
coolant velocities, and coolant-side pressures chosen. The task remained to per- 
form a final analysis of this design configuration to ensure that the desired 
features had been achieved. Before this analysis is discussed, however, a brief 
mention of the logic behind selection of water bulk temperature is presented. A 
study was conducted employing Eq. 3-4 (Table 3-9) to estimate the cooling effec- 
tiveness of water. The following values were assumed: 


200 ft/sec 


T WC = 400 F 


channel ~ 

ro ^0.05 
(Fey) 


0.06 inch 
s const ~ 2. 21 


C- = 1 Btu/lbm-R 

P 


With these values, Eq. 3-4 becomes: 


0.1961 


(Pr) 


0.6 


9 

The cooling effectiveness, Q/A, may be written; 


(3-5) 


Q/A 


h c (T WC ' T BC> 


(3-6) 


Combining Eq. 3-5 and 3-6 and introducing the value of T^ = 400 F yields: 


Q/A 


0.239 


(Pr) 


0.6 


C400 - Tbc ) 


(3-7) 


It can be seen from Eq. 3-7 that, as Tg C is increased, the effective driving po- 
tential is decreased. However, for water, the reduction of the Prandtl number 
with increased temperature overrides the loss in AT up to roughly, 200 to 250 F. 
Equation 3-7 is plotted in Fig. 3-16. It is evident that it is beneficial to heat 
the coolant water to at least 120 to 140 F. 
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A final analysis of the chamber liner was conducted assuming the coolant parameters 
shown in Fig. 3-15 and with T bC = 120 F. A two-dimensional conduction model was 
employed in place of Eq. 3-1 and 3-2. The results of that analysis are presented 
in Table 3-11. These results indicated that the first design for the chamber liner 
was more than adequate. The temperatures were quite acceptable and the cooling 
margins were greater than 50 percent. 


TABLE 3-11. FINAL CHAMBER HEAT TRANSFER RESULTS 
(T BC « 120 F) 



Throat 

Chamber 

• 2 
Q/A, Btu/in. -sec 

35 

20 

T F 

WG’ 

800 

680 

T F 

400 

430 

WC (Root) ’ 

2 

h , Btu/in. -sec-R 
g 

0.00685 

0.00414 

2 

h c> Btu/in. -sec-R 

0.090 

0.060 

V , ft/sec 

200 

153 

P„„. m , psia 

700 

538 

STATIC F 



T , F 

500 


SAT 



2 

(Q/A) b q , Btu/in. -sec 

53.7 

43.9 

Cooling Margin, 

1.53 

2.20 

(Q/A) B> 0 /CQM) 




3. 3. 2. 2 Hardware Design Description . The entire thrust chamber assembly is shown 
in Fig. 3-17, including the injector assembly previously discussed. Figure 3-17 
is the actual chamber top assembly drawing and contains all of the Rocketdyne 
drawing numbers for each part. The detail drawing of the chamber liner (presented 
as Fig. 3-18) shows all of the critical chamber dimensions (the -013 chamber was 
fabricated) . 
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The chamber liner was fabricated from NARloy (a special copper alloy) rather than 

OFHC copper to gain a high yield strength (Fig. 3-18). The chamber diameter is 

0.6 inch and the throat diameter is 0.395 inch, giving a contraction ratio 

£ = 2.307. 

c 

The liner is retained in the thrust chamber body (R1D4273, Fig. 3-17) between the 
forward closure (R1D3488) and the aft closure (R1D4272). The liner is restricted 
from moving downstream by retainer R1D4272-3 but is free to move upstream by slid- 
ing in the forward closure R1D3488. This was required to allow for thermal expan- 
sion of the lineT to avoid upset yielding. All pieces, except the liner, were 
fabricated from series 321 stainless steel. 

Cooling water is introduced through two inlets in the aft closure, R1D4272 (only 
one shown due to section) into a manifold formed between the aft closure and the 
body. The coolant enters the axial channels between the liner and the body and 
travels upstream to a manifold between the body and the forward closure (Fig. 3-17), 
The coolant exits from two outlets in the forward closure (only one shown due to 
section). In operation, the inlets and outlets were actually oriented at 90 degrees 
to one another to prevent preferential circulation. 

A pressure tap was provided in the body (R1D4273) at the downstream end of the 
cooling channels to ensure that the minimum static pressure in the cooling system 
was sufficiently high during firing to maintain the high values of T gAT required. 

Two chamber pressure taps were provided (at 180 degrees apart; only one shown due 
to section) 1.232 inches upstream of the throat. 

A simple exhaust nozzle (R1D4270) was clamped to the aft closure. The actual 
expansion contour was composed of an annular water jet that was injected around 
the hot-gas jet. This was done to reduce the cost of the chamber and was justi- 
fied due to the fact that thrust data were not to be taken with this chamber. 
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The chamber assembly is bolted together by the studs (R1C1237-019) shown in 
Fig. 3-17. With these studs in place, the chamber could be removed from the in- 
jector housing plate (R1D2964) and transported separately as an assembly. The 
chamber assembly is installed on the housing plate by studs R1C1 238-013 (not 
shown due to section), which screw into holes in the housing and protrude all the 
way through the chamber and are secured on the downstream end of the nozzle with 
nuts . 

The housing (R1D2964) was the only piece of hardware that was attached directly 
to the test stand (on plate R1D2965). The thrust chamber or the injector could 
be removed, independently, without removing the other. The injectors mount to the 
housing as an assembly so that they can be "made up" ahead of a test in the work 
shop. 

Photographs of selected components are presented in Fig. 3-19 through 3-23 and 
photographs of the body assembly (R1D4273), showing the liner in place, are pre- 
sented as Fig. 3-24 and 3-25. The entire thrust chamber assembly is shown mounted 
on the test stand (Appendix E) in Fig. 3-26. 
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5AA34- 1/7/72- SIB 


Figure 3-19. Housing (Injector End), R1D2964 


Figure 3-19. 


5AA34-1/7/72-S1A 


(Concluded) (Chamber End) 
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5AA34-1/7/72-S1F 


Figure 3-20. Forward Closure (Forward End), R1D3488 


UZ.60- 



5AA34-1/7/72-S1E 


Figure 3-21. Forward Closure (Aft End), R1D3488 
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5AA34-1/7/72-S1G 


Figure 3-22. Aft Closure (Forward End), R1D4272 
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Figure 3-23 


5AA34-1/7/72-S1H 


Aft Closure (Aft End), R1D4272 
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5AA34-2/ 3/ 72- SIC 


Figure 3-24. Body Assembly (Side View), R1D4273 
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5AA34- 2/3/72- SIB 


Figure 3-25. Body Assembly (Looking Aft) R1D4273 
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5AA24-2/22/72-S1 


Figure 3-26. Chamber Assembly on Lima Stand (R1D1243) 


3.4 EXPERIMENTAL RESULTS 


All gas/liquid experimental results are presented in this section. These include 
cold- flow mixing, cold-flow atomization, and hot-fire results. Cold- flow mixing 
results are presented in terms of the mixture ratio uniformity parameter, E m 
(Appendix F), cold- flow atomization results are presented in terms of the mass 
median droplet diameter, D, and the hot-fire results are presented in terms of 

In the section devoted to correlation of the results, cold-flow results are 


T U*- 

c 


transformed from E and D into n * and n * by means of analytical combustion 


m 


mix 


models (Appendix F) . 


vap 


The format of the presentation will be the same for all three sections and is 
listed below. (The discussion of the hot-fire results does not include a "Generali- 
zation of the Results" section.) 


1. Nominal Conditions-Baseline Comparisons 

2. Variations About AR = 3 

3. Generalization of the Results 


This format has been adopted to focus attention upon the primary objective of the 
the study, the determination of the effect of element shape upon injector perfor- 
mance at a fixed set of nominal conditons. Broad discussion of methods of gener- 
alizing the results is incoprorated in separate material at the end of each section 
and may either be studied or ignored and will not influence the general results of 
the program. These generalizations are highly interesting and important and war- 
rant incorporation in the report, however; it is suggested that they be studied 
upon a second reading of the report, if desired. The nominal, or baseline condi- 
tions are listed in Table 3-4 and the element comparison scheme has been outlined 
in Section 3.2.2 (Fig. 3-4). 
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5.4.1 Mixing Studies 


The objective of the mixing studies was to determine the effect of element geometry 
and, to some extent, operational parameters upon the mixture ratio uniformity 
parameter, E m (Appendix F) . The method for obtaining the experimental data is out- 
lined in Appendix A. Mass fluxes and mixture ratios at discrete points in the flow- 
field produced by a given element were measured by means of a two-phase sampling 
probe that was positioned 2 inches downstream of the injector element exhaust plane. 
None of the Phase II mixing results are included in this presentation, as these 
data were obtained at a sampling distance of 5 inches. Since mixing efficiency 
varies with distance, the 2- and 5-inch results would not be expected to be 
comparable 

All of the Phase III mixing cold-flow data are presented in Table 3-12. 

3. 4. 1.1 Nominal Conditions--Baseline Comparison . The mixing results for four 
baseline injector elements are presented in Fig. 3-27. Data, E m> are plotted with 
respect to actual center-post .recess in Fig. 3-27a and with respect to relative 
recess, R/y^ in Fig. 3-27b. The relative recess appears to offer the better 
representation of the results. A cross plot of Fig. 3-27b is presented as Fig. 3-28, 
wherein E is shown as a function of aspect ratio, x/y., for constant values of 
centerpost recess. (Note: Results presented at aspect ratio unity are those for 
a circular concentric tube element, not a rectangular element of aspect ratio 
unity. ) 

The representations of Fig. 3-28 show that the quality of mixture ratio uniformity 
improves dramatically with increased aspect ratio at the nominal operating condi- 
tions. This result was so interesting that a more detailed analysis of the data 
was undertaken in an attempt to explain exactly how the change in shape, alone, 
could improve the mixing to such a degree. The results of this extended analysis 
are presented in Fig. 3-29 and 3-30. 
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TABLE 3-12 NONCIRCULAR ORIFICE PROGRAM, GAS/LIQUID MIXING RESULTS (PHASE III DATA ONLY) 



Element 

Code 

No. 



(psia) 


Recess 

(in.) 



GCR-2 
GCR-2 
GST-1 
GCR-2 
GCR-2 
GCR-2 
GCR-2 
GCR-5 
GCR-S 
GCR-3 
GCR-3 
GCR-2 
GCR-4 
GCR-3 
GCR-2 
GCR-7 
GCC-2 
GCC-2 
GCC-2 
GCC-2 
GST- 2 
GCR-6 
GCR-6 
GCR-2 
GCR-2 
GCR-3 
GCC-2 


0.432 0.370 0.0617 6.0 


0.432 

0.216 

0.216 

0.432 

0.216 

0.216 

0.432 


0.370 

0.185 

0.185 

0.370 

0.185 

0.185 

0.370 


GCC-2 


0.432 

0.216 

0.324 

0.432 

0.340 

0.247 

0.154 

0.189 

0.112 

0.278 

0.159 

0.0991 

0.162 

0.0924 

0.0577 

1.024 

0.518 

0.332 

0.239 


0.370 

0.185 

0.278 

0.370 

0.278 

0.185 

0.0926 

0.127 

0.050 

0.0238 

0.119 

0.0595 

0.139 

0.0693 

0.0346 

0.878 

0.372 

0.186 

0.0930 


0.0617 
0.0309 
0.0^09 
0.0617 
0.0309 
0. 0309 
0.0617 


0. 0j317 
0.0309 
0.0^63 
0.0617 


0.0617 

0.0396 

0.0596 

0.0596 

O.Opl 

0.0231 

0.0231 

0.146 


0.146 


6.0 

4.51 

3.00 
1.50 
2.06 
0.810 

6.01 

3.01 
1.50 

6.02 

3.00 
1.50 

6.01 
2.55 
1.27 
0.637 


0.270 

0.637 


0.637 


42.5 

100.7 
42.5 
42.5 
42.5 
21.2 
31.9 

42.5 

31.9 

21.3 

10.6 

14.6 
5.74 

27.3 

13.7 
6.83 

15.9 
7.96 
3.98 

100.8 

42.7 

21.4 

10.7 


47.0 

51.0 
51.0 
51.0 

106.0 


106.0 


0.0965 

0.007 

0.0652 

0.0870 

0.1188 

0.1054 

0 

0.0475 

0.127 

0.105 

0.105 

0 

0 

0.007 

0.004 

0.007 

0.102 

0.241 

0.164 

0 

0.2375 

0.0110 

0,0110 

0.0110 

0.007 


0.00433 




1. 790 
.055 
.440 
.278 

0 

0.798 

1.539 

1.273 

1.273 

0 

0 

0.085 

0.069 

0.044 

0.636 

1.502 

1.022 

0 

2 . 

0 . 
0.133 
0.133 
0.044 


0.00433 

0.00218 

0,00316 

0.00433 


0.007 


0.00433 

0.0102 


0.044 0.0102 


22.0 

22.0 

9.62 

22.0 



0.880 

0.441 

0.880 


22.0 

9.35 

22.0 

22.0 

22.0 

43.8 

30.2 

22.0 

29.4 

44.0 

88.0 

64.1 

163.0 

22.0 

43.8 

87.8 

22.0 

44.0 

88.0 

9.28 

22.0 

43.8 
87.6 


0.880 

0.876 

0.905 

0.880 


0.880 

0.564 

0.564 

0.564 

0.316 

0.316 

0.316 

0.900 


0.900 


82.53 

63.83 

80.40 

60.69 

76.52 
80.02 
79.60 

65.48 
81.56 

77.63 
72.82 
74.59 
75.35 
49.33 
59.68 
64.37 
48.90 

49.62 
59.76 
53.89 
84.05 
57.39 

75.64 

64.62 
77.66 
68.04 

45.65 

60.62 

79.32 

81.33 

80.19 

76.28 

44.19 

75.29 
75.58 

46.18 
71.31 

77.19 

56.52 

77.29 

84.66 

82.48 



95.38 

86.33 

94.45 

95.36 

95.19 

89.03 

96.03 
94.62 
92.77 

91.54 
93.99 
77.05 
84.66 

88.15 
76.29 
76.49 
85.08 

80.44 
97.12 

83.44 

92.40 
88.27 

94.52 
90.35 
74.95 

85.54 
96.65 
98.93 

98.10 

99.15 

74.10 
94.85 
98.23 
75.21 
93.73 

98.40 

83.52 

97.11 
99.42 
99.60 
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TABLE 3-12 (Concluded) 


Element 

Code 

No. 


GCC-2 


P 

c 

(psia) 


Recess 

(in.) 


E V 

m MIX 


0.269 

0.140 

0.129 

1.09 

0.199 

0.0697 

0.129 

0.540 

0.164 

0.0348 

0.129 

0.270 

0.333 

0.239 

0.0939 

2.55 

0.213 

0.119 

0.0939 

1.27 

0.154 

0.0597 

0.0939 

! 0.636 

0.194 

0.139 

0.0548 

2.54 

0.124 

0.0696 

0.0548 

1.27 

0.0896 

0.0348 

0.0548 

0.635 

0.164 

0.0348 

0.129 

0.270 

0.149 

0.0204 

0.129 

0.158 

0.435 

0.306 

0.129 

2.37 

0.217 

0.186 

0.0310 

6.00 

0.124 

0.0930 

0.0310 

3.00 

0.0775 

0.0465 

0.0310 

1.50 

0.139 

0. 119 

0.0199 

5.98 

0.0796 

0.0597 

0.0199 

3.00 

0.0498 

0.0299 

0.0199 

1.50 

0.0812 

0.0696 

0.0116 

6.00 

0.0464 

0.0348 

0.0116 

3.00 

0.0290 

0.0174 

0.0116 

1.50 

0.119 

0.0248 

0.0939 

0.264 

0.199 

0.0697 

0.129 

0.540 


1.50 

16.0 

350 

284.0 

1.S0 

8.00 

350 

284.0 

1.50 

4.00 

350 

284.0 

0.637 

27.4 

600 

116.0 


13.7 

600 

116.0 


6.85 

600 

116.0 

if 

15.9 

7.95 

350 

A 

121.0 

121.0 

0.637 

3.98 


121.0 

1.50 

3.98 


284.0 

1.50 

2.33 


284.0 

1.50 

35.0 

350 

284.0 

0.135 

21.3 

935 

22.0 


10.6 1 

935 

22.0 


5.31 

935 

22.0 


13.6 

600 

25.0 


6.82 

600 

25.0 


3.41 

600 

25.0 


7.95 

350 

26.0 


3.98 

350 

26.0 

0.135 

1.99 

350 

26.0 

0.637 

2.84 

600 

116.0 

1.50 

7.95 

350 

284.0 


.0102 

,0240 

.0240 

.0240 

.00216 


,00216 

,0102 

,0240 


21.9 

43.8 

87.5 


21.9 

43.8 

87.6 

22.0 

44.0 

88.0 
88.0 

150.0 

10.0 

44.0 

88.0 

176.0 

44.0 

88.0 
176.0 

44.0 

88.0 

176.0 

211.0 

44.0 


0.316 

86.40 

99.59 

0.316 

83.78 

99.83 

0.316 

89.45 

99.97 

0.552 

75.62 

96.21 

0.552 

85.04 

99.45 

0.552 

81.63 

99.62 

0.316 

77.90 

96.86 

A 

85.36 

99.47 


80.18 

99.57 


90.93 

99.96 

T 

89.50 

100.03 

0.316 

79.05 

98.09 

0.909 

49.22 

75.80 

0.909 

68.23 

92.09 

0.909 

67.56 

96.98 

0.548 

51.59 

77.26 

0.548 

71.98 

94.06 

0.548 

70.83 

93.91 

0.314 

71.00 

92.18 

0.314 

64.04 

92.43 

0.314 

67.67 

97.39 

0.552 

90.81 

99.97 

0.316 

85.11 

99.86 
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Figure 3-28. Mixture Ratio Uniformity Comparison as a Function 
of Aspect Ratio (nominal conditions) 
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In these two figures, mass fraction flux (mass fraction per unit area) mff contours 

for the liquid and the gas component are shown for each of the four elements at 

zero recess. The mass fraction flux is computed by dividing the local value of 

th 

the mass flux of the 1 component, ft. /A, (liquid or gas) by the total injected 

th ^ 

mass flow of the i component, w Ti< At a point in the flowfield where the local 
mixture ratio, w L /w G is equal to the injected mixture ratio, « LT /w GT » the values 
of the mass fraction flux for the liquid and gas components would be equal. Thus, 
for an element which could produce perfect mixing, E^ = 100 percent, the liquid 
and gas mass fraction flux contours would be, everywhere, superimposed, one upon 
the other. In Fig. 3-29 constant mff contours ^w^/A) / w^ = O.sjare presented 
for the liquid and gas components. (The value mff = 0.5 was selected arbitrarily 
and the comparison could have been made at other values. Only one contour value 
was selected to avoid confusion.) These contours are obtained by plotting the 
distance from the centerline of the element, in polar coordinates, at which the 
given value of mff is found in the sampling plane (which is normal to the element 
centerline) . Here again, if the elements had produced perfect mixing, the two 
contours would lie directly upon one another. For the circular element, the 
value of 0.5 for the mff is achieved much closer to the element centerline than 
that of the gas component. There are no points in the flowfield that are "on 
mixture ratio" for this value of mff. Physically this means that the liquid is 
"coring" or remaining in the center of the flowfield and has not been spread out 
in the radial direction by the gas. 

This is a classic problem for circular concentric tube elements operating at high 
mixture ratios (>4:1) and having relatively low gas-to-liquid velocity ratios. 

The influence of element shape is quite evident in Fig. 3-29. The "warping" of 
the flowfield causes contours to approach each other or intermingle, and thus, the 
degree of mixing is increased without the necessity of increasing the gas velocity 
or decreasing the liquid velocity. 

Exactly the same result can be noted in Fig. 3-30, wherein the value of mff for 
the liquid and gas component is plotted as a function of the radial distance from 
the element centerline along a fixed ray in the collection plane (i.e., mff = f(r). 
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Figure 3-29. Effect of Liquid Port Aspect Ratio Upon Flowfield, 
R/y, = 0 (Constant mff Contours) 




6 = constant). Here, the coring problem encountered by the circular element can 
be seen plainly. The increased aspect ratio systematically reduces the "coring." 

There is another interesting potential advantage of RCTE's that can be deduced 
from Fig. 3-29. The mff contours for the highest aspect ratio element (GCR-5) 
suggest that the flowfield can be "tailored" to provide "fuel-rich" and "oxidizer- 
rich" zones that can be employed to ensure thrust chamber/injector compatibility 
in the outer zone of the chamber wall. 

Centerpost recess produces approximately the same change in the character of the 
contour plots as does increased aspect ratio. Recess has the ability to reduce 
coring . 

3. 4. 1.2 Variations About AR = 3 . Continuing with the investigation of the effect 
of shape as outlined in Fig. 3-4, the next task was to study the influence of 
element area ratio, element scale, and deviation from uniform gas annulus gap with 
aspect ratio, x^/y^, equal to 3.0. 

The effect of a reduction in the gas annulus area (with uniform gap) is presented 
in Fig. 3-31. This comparison could not be made at the nominal flowrates due to 
the fact that the gaseous nitrogen employed as a fuel simulant was already at a 
velocity of 950 ft/sec in the baseline element. A reduction in the gas area at 
the nominal conditons would have caused the nitrogen to choke. The comparison 
was made, therefore, at one-half the nominal total flowrate. The mixture ratio 
and gas density were held at nominal values, 6.0 and 0.27 lbm/ft 3 , respectively. 

The value of E m dropped off for the baseline AR = 3 element (GCR-2) as the total 
flow was reduced. At one-half nominal flowrate, the liquid velocity was 21 ft/sec 
and the gas velocity was reduced from 935 to 468 ft/sec. At these flow conditions, 
the element with reduced gas area (GCR-4) was operating with the same liquid velocity 
but with a gas velocity of 935 ft/sec. It can be seen that this increase in gas- 
to- liquid velocity ratio produced a significant increase in E . This change in 
E m is ascribed totally to the change in gas velocity, as the density and mixture' 
ratio were held constant. 
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PERCENT 



LIQUID INJECTION VELOCITY, FT/SEC 
(ALSO TOTAL FLOWRATE X 100, APPROXIMATELY) 


Figure 3-31, Effect of Variation of Gas to Liquid Velocity 
Ratio at Constant Mixture Ratio and Constant 
Gas Density (GCR-2 and GCR-4) 
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The effect of element size is shown in Fig. 3-32. Here, again, the comparison 
was not made at the nominal flowrate. The scaled-down model, GCR-7 f — - 1 - — 
reduction of GCR-2), had liquid and gas area exactly one-half those of the large 
element. Comparison was made at nominal mixture ratio, density, and velocities, 
thus the total flow in GCR-7 was exactly one-half that in GCR-2. This meant that 
the mass fluxes (w./A.) for both elements were identical. The curve of E for 

^ ^ III 

GCR-2 is shown as a function of liquid velocity. This curve is reproduced from 
Fig. 3-31. The data point for GCR-7 is plotted on the same scale. It can be 
seen that the value of E m is unaffected by size as long as the mass flows per unit 
area are preserved along with densities, velocities, and mixture ratio. This is 
a highly important result and suggests, at least at this scale, that thrust per 
element has no effect upon mixing of an individual element. This is not to say, 
however, that improved interelement mixing cannot be improved with lower thrust 
per element and thus the mixing of the injector improved overall (i.e., more 
elements] . 

The effect of variation of the gas port aspect ratio (change of x^/y^ with gas 
annulus area and centerpost configuration held constant) is shown in Fig. 3-33. 

This change produces a nonuniform gas gap. Here, values of E m for GCR-2 and GCR-3 
are plotted as functions of centerpost recess. This result, alone is relatively 
uninteresting in that it merely says that the level of mixing dropped when the 
gas port aspect ratio was changed from 1.86 to 2.347. The more interesting 
aspects of this experiment are presented in Fig. 3-34. Here, again, mff contours 
are employed to tell the story. Gas and liquid contours for mff - 0.5 are pre- 
sented for GCR-2 and GCR-3 at their respective points of centerpost recess at which 
maximum mixing level was encountered. These points are marked by arrows in Fig. 3-33. 

The contours for element GCR-2 (Fig. 3-34) suggest that there is "too much" gas 
concentraction in the y direction and that a redistribution of some of the gas 
to the X direction should tend to cause the contour lines to more closely coin- 
cide. As a matter of fact, the design of GCR-3 was based upon analysis of these 
contours and its objective was exactly that stated. What happened is obvious. 
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Figure 3-32. Effect of Element Size at Constant Density and Mixture 
Ratio (GCR-2 and GCR-7) 
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RELATIVE RECESS, R/y ] 


Figure 3-33. Effect of Variation of Gas Port Aspect 
Ratio Upon Mixing at Nominal Conditions 
(GCR-2 and GCR-3) 
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Figure 3-34. Effect of Gas Port Aspect Ratio on Flowfield Development 




The redistribution of the gas component from y to x was greatly overdone. As a 
result, the mixing uniformity dropped; however, the character of the flowfield was 
changed. This result suggests that an optimum value of X^/y^ must exist between 
the values that were tested. The inferences of these data are presented in Fig. 3-35, 
which is highly speculative, but interesting. If the trends suggested in Fig. 3-35 
are correct, then RCTE's offer yet another degree of freedom to the injector designer 
to aid in improving mixing with element geometry. 

3.4. 1.3 Generalization of the Results . Mixing uniformity cold-flow tests conducted 
with RCTE's proved to be relatively time consuming and, therefore, costly. This 
was due to the fact that the flowfields were highly nonaxisymmetric and a great 
number of sampling points had to be incorporated to obtain accurate results. On 
the other hand, cold-flow tests with the circular concentric tube element (GCC-2) 
proved to be quite economical due to the axisymmetric nature of the flowfield. 
Therefore, it was decided that an appreciation of the influence of the flow para- 
meters upon the mixing characteristics of G/L CTE's could best be obtained with the 
CCTE . An extensive set of data was obtained with element GCC-2 covering a broad 
range of gas velocity (350 to 935 ft/sec), gas density (0.135 to 1.5 ^ ), and 

liquid velocity (2 to 100 ft/sec) . This test matrix was in keeping with the state- 
ment made in Section 3.2.1, that the most practical variables of the cold-flow 

facilities were V~, p_, and V t . The results of this parametric study are presented 
b b L 

in Fig. 3-36 through 3-39. It must be pointed out that even this extensive para- 
metric variation does not show the independent effects of velocity ratio and mixture 
ratio. This is due to the fact that the entire set of results was obtained with 
one element and, thus, the area ratio, A G /A L> was not a variable. As a result, 
the mixture ratio was a dependent variable as shown in Eq. 3-8. 

MR = K (-—b- ) C3-8j 

P G V G 

where K is a constant. It has been shown (Fig. 3-31) that a variation of velocity 
ratio (by variation of area ratio) produces a change in mixing with constant density 
and constant mixture ratio. This type of comparison is not allowed by the restric- 
tions of Eq. 3-8. 
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PERCENT 



Figure 3-36. E m as a Function of V, Along Lines of Constant p G . Circular Concentric 
Tube Element. (Limited data points shown for reference; curves are ’’best 
fit" of actual data) 
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Figure 3-37. E m as a Function of Vl Along Lines of Constant p^. Circular Concentric 
Tube Element (data points not shown) 
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Figure 3-39. as a Function of V L Along Lines of Constant V G and Mixture Ratio. 
Circular Concentric Tube Element (data points not shown) 




Figures 3-36 through 3-38 show the variation of E m with liquid velocity along 
lines of constant density with gas velocity a constant in each figure. One of 
the many possible cross plots is presented in Fig. 3-39, wherein the variation 
of E m with liquid velocity is shown along lines of constant gas velocity with 
gas density everywhere constant. Lines of constant mixture ratio have been in- 
cluded in Fig. 3-38 and 3-39 to illustrate the fact that this quantity is dependent 
upon the other variables, and also to show that although the parametric roles of 
V , p„, and V. are quite clear, the interpretation of performance variations with 
constant mixture ratio throttling (V G = const) are not straightforward. 

Figures 3-36 through 3-39 show clearly that the liquid velocity is the parameter 
which most influences the mixing performance of a concentric tube element. In 
addition, the lines of constant mixture ratio in Fig. 3-38 and 3-29 show that 
for this particular element, no combination of variables can produce high levels 
of mixing as long as the mixture ratio is maintained at high values. However, 
based upon the results shown in Fig. 3-31 for a RCTE, it can be expected that a 
change in area ratio that increased the velocity ratio V^/V^ at high mixture ratio 
will produce a higher mixing level. Great care must be taken not to associate di- 
rectly the effects of changes in velocity ratio shown by these data with a change i 
velocity ratio at constant mixture ratio. These additional data are required to 
complete the study. 

These data do, however, provide certain broad design criteria. In general, one 
can conclude that "good” mixing, with zero post recess may be achieved with high 
gas velocity, high gas density, and low liquid velocity. The most important of 
these being low liquid velocity. 

3.4.2 Atomization Studies 


Results of the atomization studies are presented with the same format that was 
employed for the mixing results. Methods for obtaining data, atomization facilitie 
descriptions, and methods for dropsize analysis may be found in the Appendixes. 

The data are presented in Table 3-13. 
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TABLE 3-13 NONCIRCULAR ORIFICE PROGRAM, ATOMIZATION DATA (PHASE II RESULTS) 


Element 
Test Code 
No. No. 



1 

GCC-1 

0.0384 

0.0285 

0.00990 

2.88 

2 


0.0714 

0.0615 

0.00990 

6.21 

4 


0.0599 

0.0450 

0.0149 

3.02 

5 


0.106 

0.0915 

0.0149 1 

6.14 

8 


0.139 

0.119 

0.0200 

5.95 

9 


0.0805 

0.0605 

0.0200 

3.02 

10 

GCC-1 

0.153 

0.133 

0.0200 

6.65 

11 

GCR-1 

0.0399 

0.0300 

0.00990 

3.03 

12 

A 

0.0679 

0.0580 

0. 00990 ! 

5.86 

13 


0.0609 

0.0460 

0.0149 

3.09 

14 


0.105 

0.0905 

0.0149 

6.07 

15 

T 

0.0778 

0.0580 

0.0198 

2.93 

16 

GCR-1 

0.136 

0.116 

0.0198 

5.86 

18 

GCC-1 

0.0405 

0.0300 

0.0105 

2.86 

19 

4 

0.0705 

0.0600 

0.0105 

5.71 

20 

T 

0.0608 

0.0450 

0.0158 

2.85 

21 


0.104 

0.0890 

0.0155 

5.74 

22 

T 

0.0808 

0.0600 

0 . 0208 

2.88 

23 

GCC-1 

0.140 

0.120 

0.0200 

6.00 

25 

GCR-1 

0.0400 

0.0300 

0.0100 

3.00 

26 

4 

0.0685 

0.0585 

0.0100 

5.85 

27 

T 

0.0590 

0.0440 

0.0150 

2.93 

28 

T 

0.105 

0.0900 

0.0150 

6.00 

29 

GCR-1 

0.0810 

0.0610 

0.0200 

3.05 


GCR-2 


GCR-2 


0.432 

0.280 

0.140 

0.279 

0.140 

0.279 

0.280 

0.140 

0.431 

0.279 

0.431 


0.370 

0.240 

0.120 

0.240 

0.120 

0.210 

0.249 

0.120 

0.370 

0.240 

0.370 


0.0616 

0.0396 

0.0198 

0.0393 

0.0197 

0.0694 

0.0307 

0.0198 

0.0610 

0.0394 

0.0615 


P G 

lbm 

ft 3 

V L 

ft 

sec 

V G 

ft 

sec 

0.055 c - 

4.28 

618 

0.0555 

9.24 

618 

0.0584 

6.76 

884 

0.0584 

13.7 

884 

0.0621 

17.9 

1110 

0.0621 

9.09 

1110 

0.0621 

20.0 

1110 

0.0564 

4.51 

720 

0.0564 

8.71 

720 

0.0603 

6.91 

1010 

0.0603 

13.6 

1010 

0.0650 

8.71 

1250 

0.0650 

17.4 

1250 

0.0558 

4.51 

652 

0.0558 

9.01 

652 

0.0590 

6.76 

928 

0.0588 

13.4 

913 

0.0627 

9.01 

1140 

0.0621 

18.0 

1110 

0.0565 

4.51 

726 

0.0565 

8.79 

726 

0.0604 

6.61 

1010 

0.0604 

13.5 

1010 

0.0652 

9.16 

1260 



0.265 
0.265 
0.263 17.9 

0.171 35.9 

0.0852 17.9 

0.265 31.4 

0.260 37.2 

0.132 17.9 

0.409 55.3 

0.521 35.9 

0.808 55.3 


c Recess 
(psia) (in.) 





(PHASE III RESULTS) 


61.0 

65.2 

67.5 

39.5 

19.8 

58.9 
66.4 

32.6 

100.0 

135.0 

208.0 


0.00116 

0.00116 

0.00122 

0.00122 

0.00130 

0.00130 

0.00130 

0.00118 

0.00118 

0.00126 

0.00126 

0.00136 

0.00136 

0.00116 

0.00116 

0.00123 

0.00123 

0.00132 

0.00130 

0.00118 

0.00118 

0.00127 

0.00127 

0.00137 


144.0 
66.9 

131.0 

64.3 

62.4 

122.0 

55.8 
160.0 

82.6 

147.0 

74.6 

144.0 

71.8 

145.0 

72.3 

137.0 

68.3 

127.0 

61.9 

161.0 

82.7 

154.0 

75.4 

137.0 



0.929 

0.929 

0.572 

0.572 

0.832 

0.832 

1.06 

1.06 

0.514 

0.514 

0.753 

0.740 

0.961 

0.929 

0.577 

0.577 

0.837 

0.837 

1.07 


P G V G 2 

(psia) 

D 

(microns) 

D/Xi 

4.57 

157 

0.0386 

4.57 

374 

0.0920 

9.84 

133 

0.0327 

9.84 

268 

0.0659 

16.7 

297 

0.0731 

16.7 

166 

0.0408 

16.7 

393 

0.0967 

6.31 

108 

0.0733 

6.31 

151 

0.102 

13.4 

107 

0.0726 

13.4 

126 

0.0855 

21.9 

120 

0.0815 

21.9 

156 

0.105 

5.11 

167 

0.0411 

5.11 

197 

0.0485 

10.9 

152 

0.0374 

10.5 

183 

0.0450 

17.8 

134 

0.0330 

16.7 

188 

0.0463 

6.43 

123 

0.0835 

6.43 

186 

0.126 

13.5 

119 

0.0808 

13.5 

152 

0.103 

22.3 

124 

0.0842 



0.00556 

17.3 

0.785 

52.6 

490 

0.234 

0.00555 

17.2 

0.488 

21.7 

520 

: 0.248 

0.00551 

17.3 

0.241 

5.48 

700 

0.334 

0.00358 

26.5 

0.775 

33.2 

420 

0.200 

0.00179 

53.3 

0.780 

16.8 

205 

0.0978 

0.00555 

34.5 

0.900 

66. 8 

210 

0.100 

0.00545 

13.1 

0.378 

13.3 

490 

0.234 

0.00276 

34.5 

0.489 

10.9 

420 

0.200 

0.00858 

11.1 

0.487 

33.4 

545 

0.260 

0.0109 

8.71 

0.241 

10.9 

615 

0.293 

0.0169 

5.68 

0.243 

17.2 

610 

0.291 
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TABLE 3-13 (Continued) 


948 

616 

296 

945 

908 

1060 

471 

979 

638 

321 

963 

971 

1090 

480 

958 
610 
301 
934 
602 
296 
922 
925 

1040 

465 

937 

938 
938 
937 
940 
945 
931 
925 
933 
975 
971 
966 
950 
944 
942 
947 
315 

959 
1010 

953 


c 

(psia) 

61.0 

65.2 

67.5 

39.5 

19.8 

58.9 

66.4 
61.0 
65.2 

67.5 

39.5 

19.8 

58.9 

66.4 
61.0 
65.2 

67.5 
61.0 
65.2 

67.5 

39.5 

19.8 

58.9 
66.4 
61.0 


Recess 

(in.) 


Mach 

No. 


(psia) I (microns) 


0 

0.189 

0.189 

0.189 

0 


0 

0.0324 

0.0324 

0.0625 

0.0625 

0.0965 

0.0965 

0.126 

0.126 

0.0233 

0.0552 

0.0552 

0.0844 

0 


0 

0.105 

0.105 

0.105 


0 

0.393 

0.393 

0.758 

0.758 

1.17 

1.17 

1.53 

1.53 

0.392 

0.928 

0.928 

1.42 

0 


0.00572 

0.00543 

0.00576 

0.00355 

0.00189 

0.00551 

0.00553 

0.00557 

0.00549 

0.00540 

0.00363 

0.00178 

0.00560 

0.00567 

0.00547 

0.00561 

0.00574 

0.00559 

0.00559 

0.00574 

0.00366 

0.00184 

0.00567 

0.00567 

0.00563 

0.00563 

0.00563 

0.00561 

0.00560 

0.00562 

0.00565 

0.00569 

0.00555 

0.00525 

0.00525 

0.00531 

0.00540 

0.00529 

0.00530 

0.00526 

0.00518 

0.00525 

0.00319) 

0.001731 


17.6 

0.787 

52.9 

505 

0.334 

17.7 

0.482 

21.2 

460 

0.304 

17.0 

0.235 

5.21 

540 

0.357 

27.1 

' 0.768 

32.6 

345 

0.228 

52.1 

! 0.762 

16.1 

178 

0.117 

35.0 

0.885 

64.7 

202 

0.134 

13.0 

0.369 

12.6 

655 

0.433 

17.7 

I 0.803 

55.0 

400 

0. 135 

17.8 

0.502 

23.0 

460 

0.155 

17.9 

0.246 

5.73 

478 

0.162 

26.8 

0.791 

34.7 

395 

0.133 

54.2 

0.791 

17.3 

190 

0.0642 

34.8 

0.913 

68.7 

200 

0.0676 

12.9 

0.380 

13.5 

498 

0.168 

17.3 

0.778 

51.7 

284 

0.136 

17.0 

0.485 

21.5 

300 

0.143 

16.8 

0.238 

5.34 

313 

0.149 

7.36 

0.767 

50.3 

425 

| 0.203 

7.31 

0.478 

20.9 

415 

i 0.198 

7.20 

0.234 

5.20 

410 

0.196 

11.1 

0.761 

32.0 

385 

0. 184 

22.5 

0.765 

16.2 

280 

0.134 

14.5 

0.877 

63.5 

335 

0.160 

5.45 

0.369 

12.6 

270 

0.129 

16.9 

0.772 

50.9 

425 

0.203 

17.0 

0.773 

51.0 

415 

0.198 

17.0 

0.772 

50.9 

410 

0.196 

16.9 

0.770 

50.6 

385 

0.184 

17.0 

0.772 

50.9 

280 

0.134 

17.1 

0.777 

51.6 

335 

0.160 

16.8 

0.768 

50.4 

270 

0.129 

16.7 

0.766 

50.1 

333 

0.159 

17.4 

0.764 

49.8 

375 

0.248 

18.1 

0.775 

51.3 

325 

0.215 

18.1 

0.772 

51.0 

330 

0.218 

18.0 ! 

0.772 

51.0 

273 

0.181 

17.7 

0.766 

50.2 

350 i 

0.232 

54.2 

0.745 

48.5 

145 

0.0959 

108.0 

0.745 

48.4 

150 

0.0993 

72.5 

0.746 

48.6 

140 

0.0926 

17.6 

0.236 

5.29 

609 

0.291 

17.3 

0.763 

49.7 

311 

0.148 

28.3 

0.781 

33.7 

212 

0.101 

53.1 

0.764 

16.1 

175- 

0.0835 
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TABLE 3-13 (Concluded) 


Test 

No. 

Element 

Code 

No. 

*T 

lbm 

sec 

Hj 

W G 

lbm 

sec 

*L 

*G 

P L 

lbm 

ft 3 

P G 

lbm 

ft 3 

i 

V G 

ft 

sec 

P 

c 

(psia) 

Recess 

(in.) 

1 


KB 

V G 

V L 













56 

GCR-2 

0.277 

0.210 

0.0668 

3.14 

47 

.7 

0.250 

31.4 

1080 

58.9 



1.27 

0.00524 

34.7 

57 

GCR-2 

0.279 

0.249 

0.0299 

8.33 

J 


0.252 

27 

.2 

482 

66.4 

1 

1 

1.27 

0.00528 

13.0 

58 

GCC-2 

0.432 

0.370 

0.0616 

6.01 



0.267 

55 

.3 

938 

61 

.0 

0.00751 


0.0468 

0.00560 

17.0 

59 

GCC-2 

0.431 


i 

0.0610 

6.07 



0.260 



954 


i 

0.0540 


0.337 

0.00545 

17.3 

60 

GCC-2 

0.431 



0.0614 

6.02 



0.271 



923 



0.101 


0.632 

0.00568 

16.7 

61 

GCC-2 

0.431 



0.0610 

6.07 



0.262 



947 



0.156 


0.971 

0.00549 

17.1 

62 

GCR-3 

0.432 



0.0619 

5.98 



0.274 



933 



0.00231 


0.0279 

0.00574 

16.9 

63 

GCR-3 

0.432 



0.0619 

5.98 



0.271 



944 



0.0488 


0.592 

0.00567 

17.1 

64 

GCR-3 

0.432 



0.0619 

5.98 



0.271 



944 



0.101 


1.22 

0.00567 

17.1 

65 

GCR-3 

0.432 

0.370 

0.0617 

6.00 



0.269 



947 



0.151 


1.83 

0.00563 

17.1 

66 

GCR-7 

0.216 

0 . 

185 

0.0312 

5.93 



0.268 



941 



0.0030C 


0.0517 

0.00562 

17.0 

67 

GCR-7 

0.216 

0.185 

0.0313 

,5.92 



0.272 



931 



0.0353 


0.609 

0.00569 

16.8 

68 

GCR-7 

0.216 

0 . 

185 

0.0312 

5.94 



0.264 



953 



0.0645 


1.11 

0.00554 

17.2 

69 

GCR-7 

0.216 

0 . 

L85 

0.0312 

5.93 



0.268 



943 



0.104 


1.81 

0.00561 

17.0 

70 

GCR-6 

0.432 

0.370 

0.0616 

6.01 



0.271 



964 



0.0550 


0.567 

0.00567 

17.4 

71 

GCR-6 

0.431 

0.370 

0.0610 

6.07 



0.266 

55 

.3 

972 

61.0 

0.0550 


0.472 

0.00557 

17.6 

72 

GST-2 

0.140 

0 . 

120 

0.0197 

6.10 


f 

0.259 

41 

.2 

311 

67.5 

0 


0 

0.00543 

7.55 

73 

GST-2 

0.140 

0 . 

120 

0.0196 

6.11 

47 

.7 

0.0845 

41 

.2 

950 

20.0 

0 


0 

0.00177 

23.1 



0 v — 

Mach P G G D 

No. (psia) (microns) 


0.880 

0.369 


171 

555 

0.0816 

0.265 

0.771 

50.7 

381 

0.0935 

0.773 

51.1 

372 

0.0913 

0.764 

49.8 

320 

0.0785 

0.771 

50.7 

332 

0.0815 

0.776 

51.4 

326 

0.156 

0.780 

52.0 

402 

0.192 

0.780 

52.0 

330 

0.157 

0.780 

52.0 

392 

0.187 

0.775 

51.3 

288 

0.195 

0.771 

50.8 

270 

0.183 

0.779 

51.9 

284 

0.193 

0.776 

51.4 

182 

0.123 

0.797 

54.3 

370 

0.125 

0.796 

54.2 

338 

0.114 

0.239 

5.40 

400 

0.525 

0.767 

16.5 

177 

0.232 
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Atomization results are interpreted in terms of the mass median droplet diameter, 

D, obtained from a given sample of "frozen" wax droplets. The value of T) is merely 
a shorthand method of evaluating the atomization characteristics of an injector 
element. Unfortunately, D by itself does not completely describe the nature of a 
given droplet sample. The droplet size distribution function is required in addi- 
tion to D to perform a complete vaporization limited performance calculation (see 
Appendix F) , A typical droplet diameter distribution is shown in Fig. 3-40, wherein 
the cumulative weight fraction of the wax droplet sample is shown as a function of 
the normalized droplet diameter, D/D. By definition, these distribution functions 
pass through the value of 0.5 for cumulative weight fraction when D/D = 1. The 
data shown in Fig. 3-40 represent the results of three different atomization experi- 
ments conducted with two different RCTE's (see tests 5, 24, and 54 in Table 3-13). 
One distribution function fits the data for these three tests quite well. This 
distribution function was employed for all performance analysis throughout the 
remainder of the study. 

Phase II atomization data have been incorporated into this discussion. 

3. 4. 2,1 Nominal Conditions — Baseline Comparison. The nominal conditions--baseline 
comparison is shown in Fig. 3-41. Here the droplet mass median diameters normalized 
with respect to the characteristic size of the element are plotted as a function of 
the relative recess, R/y^, of the liquid centerpost. Two basic conclusions can be 
drawn from these results: (1) droplet diameter is reduced with increased center- 

post recess, and (2) relative droplet diameter increases with increased aspect 
ratio. The effect of recess is well known and has been documented many times in 
the past. The effect of aspect ratio, however, is quite surprising and warrants 
further discussion. 

A cross plot of these data is presented in Fig. 3-42, wherein the relative dropsize 
is shown as a function of aspect ratio along lines of constant relative recess. 

Two important items concerning Fig. 3-42 must be kept in mind when evaluating the 
data: (1) relative dropsize is shown, not dropsize, and the value of y. is con- 

tinually decreasing with increased aspect ratio; and (2) although these curves 
are drawn to the aspect ratio unity, the data at unity were obtained with a circular 


R-9270 

87 



R-9270 



Figure 3-40. Typical Dropsize 
Concentric Tube 



Distribution Function. Gas/Liquid Rectangular 
lement (GCR-2), and Circular Concentric 






Figure 3-41. Normalized Droplet Diameter as a Function of 
Center-Post Relative Recess for the Four 
Baseline Elements at Nominal Conditions 
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Figure 3-42. Normalized Droplet Diameter as a Function of 
Aspect Ratio for Constant Values of Relative 
Recess (baseline elements at nominal conditions) 
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element, not a RCTE with aspect ratio unity, and should not be expected to fall 
on the curves . 


It was found that a relatively simple algebraic function could be employed to 
correlate the RCTE dropsize results. This function is presented as Eq. 3-9: 



where K = f (recess). 


(3-9) 


The curves that are drawn through the data in Fig. 3-42 were generated with 
Eq. 3-9. The implication that these curves fall below the CCTE data at aspect 
ratio unity is purely speculation and has not been shown. It is interesting to 
note that Eq. 3-9 states that, geometrically, dropsize is a function of the 
term: 


(Note: 



(3-10) 


It must be remembered that the liquid port area, X^, was held constant 
as aspect ratio was varied.) 


It is interesting to speculate upon what a correlation of droplet diameter based 
upon the hydraulic diameter of the liquid orifice would yield as a functional 
relationship. If it had been assumed that the correlating parameter was D/D 

H 

where = hydraulic diameter, the following would have been the result: 


D = 


f 


*1 


AR 

1 + AR 


(3-11) 


The similarity between this function and Eq. 3-10 is striking. However, hydraulic 
diameter alone is not sufficient to describe the data. The difference between 
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Eq. 3-10 and 3-11 also point out the danger of adopting complex combinations of 
variables as correlating parameters early in an experimental study. 


The variation of actual dropsize with aspect ratio is shown as a function of 
aspect ratio in Fig. 3-43. The curve through these data was obtained from Eq. 3-10 
by adding the additional restriction that the product X^y is a constant and is 
equal to the liquid area, Aj. This restriction adds the additional relation: 



Combination of Eq. 3-9 and 3-11 yields: 



where 


(3-12) 


(3-13) 


K = 0.310 



513 microns 


AR = X 1 /y 1 


These results can be interpreted quite simply from the implications of Eq. 3-10. 
This relation states that the characteristic dimension for atomization is a com- 
plex quantity at low aspect ratios and approaches y^ as AR -*■ 00 . Physically, this 
is reasonable. At low aspect ratios, although y^ is reduced by an increase in 
aspect ratio, X^ is increased. Evidently, the increase in the X dimension is more 
detrimental to atomization than the decrease in y^ is beneficial for atomization. 
As the aspect ratio increases to larger and larger values, the characteristic 
thickness of the liquid jet quite naturally approaches the value of y^. 


R-9270 

92 



MICRONS 


1000 



Figure 3-43. Variation of Wax Droplet Diameter With Aspect 
Ratio (baseline elements, nominal conditions) 
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3. 4. 2. 2 Variations About AR = 3. Only two of the three variations about the 
aspect ratio three baseline element were investigated during the atomization 
study (element scale and nonuniform gas gap) . (Atomization experiments with 
element GCR-4, one-half baseline gas area, were not conducted due to exhaustion 
of the budget allocated to cold-flow atomization testing before completion of 
this element). 

The effects of variation of element size and gas port aspect ratio are presented 
in Fig. 3-44. The effect of size is shown in Fig. 3-44a, wherein the relative 
dropsize is presented as a function of actual recess for both GCR-2 (the larger 
element) and GCR-7 (the smaller element). What is shown in this figure is the 
reduction in dropsize produced by a / 2 reduction in the characteristic size 
of the element. If the characteristic dimension, y^, had been sufficient to des- 
cribe the droplet variation the points for GCR-7 would have fallen directly on 
the curve for GCR-2. At zero recess the GCR-7 element actually provided a drop- 
size reduction of roughly 1.2 /T . 

This result suggests that the effect of element size upon dropsize is stronger 
than the first power of a characteristic dimension. However, the limited amount 
of data plus the fairly large degree of uncertainty associated with dropsize data 
preclude any valid speculation as to the second-order effects of element size. 

The effect of the gas port aspect ratio is shown in Fig. 3-44b. At zero recess, 
the element having a nonuniform gas gap (GCR-3) produced a D significantly smaller 
than the baseline element. However, the dropsizes for the two elements approach 
one another at moderate recess. There appears to be considerable scatter in these 
data such that a firm conclusion as to the gas gap effect cannot be drawn. The 
results seem quite rational except for the data point at the greatest recess for 
element GCR-3. 

Atomization experiements with element GCR-4 (1.2 nominal gas area) were not con- 
ducted as discussed earlier. However, it can be postulated that an improvement 
in atomization would have been realized with this element at the nominal flowrate 
and mixture ratio. This improvement would be attributable to the doubling of the 
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gas velocity at constant liquid velocity and gas density. The magnitude of the 
improvement, however, is unknown for these elements (see Ref. 7 for other elements) . 

3. 4. 2. 3 Generalization of the Results . The bulk of the parametric atomization 
data was obtained with the baseline AR = 3 element (GCR-2) . As a result, these 
data were employed in an attempt to generalize the results to operating conditions 
other than the nominal set. The influence of a parametric variation of gas velocity, 
gas density, and liquid velocity at mixture ratio equal to 6.0, is depicted in 
Fig. 3-45, In this figure, the relative dropsize, D/ y^, is shown as a function of 
liquid injection velocity along lines of constant density and constant gas velocity. 
The parametric variations of V^, V^,, and p G are not independent and are subject to 
the restrictions of Eq. 3-8 with mixture ratio a constant. It was found through 
further analysis of data at other mixture ratios that the effect of mixture ratio 
could be incorporated by casting the results shown in Fig. 3-45 into the form shown 
in Fig. 3-46. Here, D/y^^ is plotted as a function of p G V G along lines of constant 
liquid velocity. 

The restrictions of Eq. 3-8 still apply here; however, the mixture ratio is allowed 

to vary in addition to the other parameters. A physical interpretation of the data 

in this form is quite clear. With a fixed liquid velocity, the droplet size is 

2 

significantly reduced by increasing P G V G . However, as the liquid velocity is 
increased, the value of p G V ^ required to produce the same droplet diameters is 
increased. The fact that the product p G improves the correlation of the data 
suggests strongly that the Weber number is a key parameter for the description of 
the atomization process (see Table 3-3) . 

This method of presentation is extended to the other baseline elements in Fig. 3-47 
through 3-49. In Fig. 3-47, the atomization results for both the Phase II (GCC-1) 
and the Phase III (GCC-2) CCTE's are presented. This figure is significant in that 
the direct correlation between low density, "open-air" atomization and high-density 
atomization can be accomplished. These data (Table 3-13) incorporate large varia- 
tions in gas density (0.055 to 0.27 lb/ ft 3 ), liquid velocity (4 to 55 ft/sec), and 
gas velocity (300 to 1000 ft/sec). 
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Figure 3-45. Effect of V L , V G> and p G Upon Dropsize at Constant 
Mixture Ratio for the Baseline AR=3 RCTE (GCR-2) 
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Figure 3-47. Variation of Dropsize With and Vl. Baseline 

Circular Element. (Note: Solid lines are estimates 

based upon the limited data shown) 
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A very interesting result is shown in Fig. 3-49 in which the dropsize data from 
both the Phase II (GCR-1) and the Phase III (GCR-5) AR = 6 elements are presented. 
Here again, the Phase II and III data correlate well except for a set of points 
which has been pointed out in the figure. It is believed that these points con- 
stitute limiting dropsizes. That is to say that a dropsize can be achieved, 
regardless of how low the liquid velocity is or how high the value of p G V G is, 
such that a smaller droplet cannot be realized. The smallest dropsizes shown 
in Fig, 3-49 are on the order of 100 to 110 microns. It is reasonable to suppose 
that such a limit does exist. As the droplets become smaller and smaller, they 
are more easily accelerated by the gas stream (the acceleration is proportional 
to 1/D) . If the droplets are accelerated rapidly to the gas velocity, the A V 
between the droplets ar.d the gas goes rapidly to zero and further droplet breakup 
is halted. This would not be the case, however, if the droplet were held in place 
until breakup occurred. In that case, smaller droplets could be obtained from 
droplets less than 100q with increased gas velocity. 

Further dropsize comparisons for the baseline elements are presented in Fig. 3-50 

2 

through 3-52 as functions of p G V G for three different liquid velocities. 

Finally, the effect of liquid velocity variation upon dropsizes with centerpost 
recess is shown in Fig. 3-53. Once again, the limiting dropsize phenomenon appears 
in these data. Actually, the dropsize at the breakoff point is approximately 200y. 
This is a rather large dropsize and smaller droplets should have been achieved for 
this element. The premature breakoff at this point has not been explained. 

3.4.3 Multi-Showerhead Triplet Results 

The basic multi-showerhead triplet (MST) element mixing and atomization results 
are discussed briefly in this section (for data see Tables 3-12 and 3-13). The 
specification for the two MST elements tested (GST-1 and GST-2) is presented in 
Table 3-5. A discussion of the geometry of these elements may be found in 
Section 3.3.1. 
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Figure 3-50. Baseline Element Dropsize Comparison 
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Figure 3.53. Dropsize Variation With V L Showing the Effect of 
Center-Post Recess (GCR-2) 
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A face pattern view of the two MST elements is shown in Fig. 3-54. In each case, 
two rectangular gas jets are impinged at an included angle of 60 degrees upon a 
set of showerhead liquid orifices. The total gas and liquid areas are equal, 
respectively, for the two elements. 

The objective of the test series was to determine the effect of the size and number 
of the liquid jets upon mixing and atomization. The basic results are presented 
in Fig. 3-55. It can be seen that the number of orifices had little effect on the 
mixing characteristics of the element. However, dropsize was greatly improved 
with the larger number of orifices (smaller jets) at low gas velocity. At high 
gas velocity the jet size had no effect. 

The truly interesting aspect of these data is that the mixing uniformity is rela- 
tively high and the droplet diameters are quite small. Further atomication charac- 
teristics of the 28-oridice element (GST-1) are presented in Fig. 3-56. Dropsize 
is shown as a function of liquid velocity along lines of constant gas velocity 
and gas density. The trend with constant gas density suggests that there is an 
optimum operating point for atomization. This is logical in view of the fact that 
this element is subject to agglomeration of the liquid jets if the gas velocity 
is not high enough. The diameter of each jet is already 500 microns. This is as 
small as some of the droplets produced by the concentric tube elements. Evidently, 
below the optimum operating point the benefit of reduced liquid velocity is over- 
come by the agglomeration produced by the ’'pushing together" of the liquid jet by 
a low gas velocity. 

In all, these are very interesting elements and should be investigated in more 
detail. The brief mention of their results in this report stems not from a lack 
of technical interest but rather from limitations dictated by the scope of the 
program. 

These elements were not included in the hot-fire experiments under this study. 
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Figure 3-55. Atomization and Mixing Results for the Multi-Showerhead 
Triplet Elements (GST-1 and GST-2) 
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Figure 3-56. Variation of Dropsize for the GST-1 Element as a 
Function of Liquid Velocity Along Lines of 
Constant V„ and p_ 
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3.4.4 Hot-Fire Studies 


Results of the hot-fire studies are presented in Table 3-14. Two characteristic 
velocity efficiency values are presented ( r| c * RAW > was obtained directly from the 
test data with the only correction being an adjustment of static chamber pressure 
to throat stagnation pressure under the assumption of isentropic flow from the 
static pressure tap to the throat). For the contraction ratio 2.3:1 chamber 
employed, this correction was: 

P = P x (1.0418) (3-14) 

TS c STATIC 


The value of C* efficiency labeled n * 

c corr 

according to the following relationship: 


has been corrected for bulk heat loss 


where 


where 


= n. 


X f, 


corr 


RAW 


HL 



AQ = total heat loss, Btu/sec 

» 

= total chamber flowrate 

C p = chamber gas specific heat, Btu/lbm-R (shifting) 
*s 

T c = theoretical combustion flame temperature, R . 


(3-15) 


(3-16) 


In general, a bulk heat loss correction is inappropriate for most large rocket 
engine applications. However, in this case, the chamber diameter is extremely 
small and the chamber length- to-diameter ratio is quite large (L/D ~ 6.7). There- 
fore, it was assumed that the bulk heat loss was a valid, first-order correction. 
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TABljE 3-14 NONCIRCULAR ORIFICE PROGRAM GAS/LIQUID HOT-pIRE RESULTS (PHASE IV) 


Element T 
Test Code lbm 
No. No. 



*G 

MR 

P 

w. 


lbm 

L 

i 

* G i 


sec 

f 



G c Oxid 

ft (psia) Recess Inj . 
sec Static (in.) (°F) 


T T 
Oxid Fuel 1 
Inj. Inj.) 
(°F) (OF); 



Mach 

No. 


GCC-2 


GCC-2 

GCR-8 


GCR-8 

GCR-4 


GCR-4 

GCR-7 


GCR-7 

GCR-10 

GCR-10 

GCR-10 

GCR-9 

GCR-9 

GCR-9 

GCR-9 


0.350 

0.463 

0.397 

0.488 

0.539 

0.587 

0.370 

0.351 

0.382 

0.503 

0.363 

0.364 

0.224 

0.494 

0.216 

0.248 

0.377 

0.209 

0.351 

0.449 

0.399 

0.307 

0.276 

0.400 

0.491 

0.273 

0.368 

0.500 

0.439 

0.304 

0.400 

0.271 

0.327 

0.290 

0.406 

0.330 

0.309 

0.411 

0.550 

0.490 



0.340 

0.419 

0.469 

0.508 

0.321 

0.301 

0.332 

0.434 

0.278 

0.301 

0.195 

0.424 

0.185 

0.197 

0.327 

0.179 

0.301 

0.385 

0.331 

0.268 

0.237 

0.345 

0.424 

0.233 

0.317 

0.431 

0.378 

0.262 

0.346 

0.228 

0.282 

0.238 

0.354 

0.288 

0.267 

0.356 

0.480 

0.426 


0.0579 

0.0574 

0.0569 

0.0689 

0.0700 

0.0785 

0.0490 

0.0496 

0.0495 

0.0693 

0.0853 

0.0626 

0.0298 

0.0699 

0.0305 

0.0505 

0.0505 

0.0304 

0.0501 

0.0639 

0.0688 

0.0396 

0.0396 

0.0553 

0.0668 

0.0399 

0.0509 

0.0693 

0.0606 

0.0419 

0.0539 

0.0436 

0.0449 

0.0521 

0.0524 

0.0425 

0.0422 

0.0550 

0.0704 

0.0633 



5.98 

6.08 

6.70 
6.48 
6.55 
6.07;, 

6.71 
6.26 

3.26 
4.81' 
6.S4 
6.08 
6.08 
3.91 
6.47 
5.88 
6.01 
6.03 
4.80| 
6.75 
5.98 
6. 23 I 
6.361 
5.85 

6.23 

6.22 

6.24 

6.26 
6.42 
5.22 
6.29 
4.57 
6.75 
6.80 
6.33 
6.47 
6.81 
6.74 


69.5 

73.7 
73.3 
73.3 

73.2 

72.8 

69.3 

71.4 

71.6 
72.0 

71.5 

72.0 

70.0 

71.9 

69.0 

68.2 

70.5 

69.3 

70.3 

71.2 

71.4 

71.1 
71.1 

70.8 

71.1 

71.5 

71.9 

71.6 

71.3 

71.6 

71.8 

63.3 

65.7 

63.9 

63.2 

58.9 
60.6 

60.4 

63.8 
63.6 


0.191 

0.231 

0.206 

0.250 

0.271 

0.300 

0.214 

0.205 

0.219 

0.290 

0.239 

0.225 

0.128 

0.293 

0.132 

0.165 

0.215 

0.133 

0.220 

0.285 

0.263 

0.180 

0.167 

0.244 

0.298 

0.172 

0.226 

0.300 

0.267 

0.189 

0.245 

6.188 

0.219 

0.183 

0.225 

0.183 

0.173 

0.229 

0.298 

0.265 


30.0 
39.3 

33.1 
40.8 

45.7 

49.8 

33.2 

30.2 

33.2 

43.2 

27.9 

30.0 

19.9 

42.3 

19.2 

20.7 

33.2 

18.5 

30.7 

38.8 

33.2 

26.8 

23.7 

34.7 

42.6 

46.5 
62.9 

85.8 

75.7 

52.3 

68.7 

51.3 

61.3 

26.7 

40.1 

34.0 

31.6 

42.3 

53.9 

48.1 


534.7 

651.8 

583.0 

707.5 

769.4 

851.9 
617.8 

594.1 

635.7 

839.7 

682.6 

650.7 

372.5 

841.0 

381.1 

472.4 

618.1 

385.0 

635.8 
821.3 

756.9 

516.0 
473.8 

694.7 

847.0 

486.0 

637.8 

844.2 

751.7 

526.1 

683.7 

528.0 

616.1 

516.2 

638.7 

513.7 

497.5 

659.0 

860.5 

764.0 


0 

0.0779 

0.0779 

0.0779 

0.1164 

0.1164 

0.1164 

0.1164 

0.1641 

0 


0 

0.0536 

0.0536 

0.104 

0.104 

0 

0 

0.0123 

0.0110 

0.0110 

0.0110 

0.0110 


284.2 

309.0 

307.0 
305.8 

304.7 

301.7 

282.5 

295.1 

296.1 

296.8 

295.2 

298.2 
288.0 

295.9 

282.7 

277.3 

289.4 

284.2 

288.2 

292.0 

293.6 

293.8 

293.9 

290.6 

291.0 

296.4 

297.5 

294.3 
292.8 

296.4 

296.7 
-250.2 
-261.1 
-254.4 
-250.0 
-232.3 
-239.5 
-236.9 
-251.2 
-250.8 


0.9442 0 
0.9442 0 
0.9442 0 
1.411 0 

1.411 0 

1.411 0 

1.411 0 

1.989 0 

0 0 


0.9241 0 
0.9241 0 
1.793 0 

1.793 0 

0 0 

0 0 

0.2067 0 
0.0944 0 
0.0944 0 
0.0944 0 
0.0944 0 


00275 41. 
00313 25. 
00281 33. 
00341 27. 
00370 23. 
00412 21. 
00309 28. 
00287 32. 
00306 27. 
00403 22. 
00334 52. 
00313 37. 
00183 47. 
00408 23. 

00191 48. 
00242 60. 
00305 28. 

00192 50. 
00313 30. 
00400 23. 
00368 32. 
00253 66. 
00235 80. 
00345 52. 
00419 42. 
. 00241 40. 
,00314 29. 
,00419 21. 

00374 24. 
,00264 34. 
,00341 25. 
.00297 36. 
,00333 27. 
,00286 44. 
.00356 24. 
,00311 27. 
.00285 32. 
.00379 24. 
.00467 18. 
.00417 21. 


2 0.290 

8 0.237 

9 0.262 

5 0.262 

0 0.245 

3 0.248 

1 0.216 

6 0.227 

7 0.212 
6 0.225 

2 0.338 

8 0.261 

3 0.218 

0 0.225 
9 0.217 

1 0.289 
9 0.222 

2 0.214 

2 0.214 
6 0.211 
1 0.247 

3 0.413 
9 0.447 
8 0.427 
6 0.422 
3 0.438 

0 0.426 

8 0.438 

3 0.430 

4 0.423 

9 0.418 

6 0.440 

1 0.389 
1 0.275 
1 0.225 
6 0.220 

7 0.239 
1 0.236 
6 0.231 
0 0.234 


.7.90 

8.23 
7.07 
8.26 
9.09 
9.53 

10.7 
10.5 
11.0 

12.8 
11.2 
11.0 

7.66 

13.4 
10.9 

14.5 

11.2 

12.4 

17.8 

21.6 

22.7 

12.6 

12.3 

14.9 

16.7 

11.0 

13.0 

13.8 
13.8 

12.5 

14.8 

18.7 

19.1 

12.2 

13.7 

9.24 
9.05 

11.3 

13.0 

12.1 


76.1 1. 

75.4 1. 

75.6 1. 

74.8 1. 

75.4 1. 

75.9 1. 

87.8 1. 

87.5 1. 

88.1 1 . 

86.5 1. 
89.4 1. 

88.3 1. 

87.7 1. 

87.7 1. 

91.7 1. 

91.9 1. 

85.9 1. 

94.7 1. 

93.3 1. 

94.0 1. 

92.0 1. 

87.4 1. 

86.7 1. 
88.2 1. 

87.9 1. 

89.5 1. 

88.1 1 . 

85.5 1. 
87.0 1. 
88.2 1. 

87.6 1. 
95.4 1. 

95.9 1. 

85.3 1. 

81.4 1. 
81.2 1. 

82.5 1. 
82.4 1. 
81.2 1. 

80.9 1. 
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Test 

Element 

Code 

No. 

No. 


60 

GCC-2 

61 



62 

63 

65 

66 



67 

GCC-2 

68 

GCR-10 

69 

GCR-10 

70 

GCR-10 

71 

GCR-3 

72 

GCR-3 

73 

GCR-3 



0.362 

0.313 

0.0492 

0.446 

0.386 

0.0594 

0.530 

0.459 

0.0707 

0.524 

0.453 

0.0709 

0.443 

0.383 

0.0599 

0.517 

0.448 

0.0688 

0.532 

0.461 

0.0704 

0.432 

0.376 

0.0564 

0.555 

0.477 

0.0783 

0.465 

0.404 

0.0619 

0.447 

0.389 

0.0585 

0.503 

0.434 

0.0688 

0.573 

0.497 

0.0751 


FOLDOUT frams 


TABLE 3-14 (Concluded) 


P G 

v. 

V r 

p 


lbm 

L 

b 

c 


V 

ft 

ft 

(psia) 

Recess 

ft 3 

sec 

sec 

Static 

(in.) 


T 

Fuel 

R 

Inj • 
(°F) 

Bza 


V 

RAW 

f 

n * 

c 

(%) 

f HL 

CORR. 


0.214 

0.255 

0.299 

0.298 

0.264 

0.298 

0.289 

0.235 

0.299 

0.251 

0.222 

0.254 

0.272 


934 

619.0 

0.153 

-249.8 

82.6 

0.9539 

0.00339 

948 

736.4 

0.153 

-261.7 

83.2 

0.9539 

0.00389 

963 

861.1 

0.153 

-258.9 

82.5 

0.9539 

0.00458 

969 

858.1 

0.153 

-264.1 

82.0 

0.9539 

0.0044S 

926 

753.3 

0.2156 

-267.3 

77.3 

1.344 

0.00396 

942 

852.2 

0.2156 

-269.8 

78.7 

1.344 

0.00443' 

993 

834.1 

0.0905 

-269.5 

84.1 

0.5642 

0.0043C 

968 

683.5 

0.0123 

-249.2 

86.9 

0.2067 

0.00372 

1057 

889.6 

0.0123 

-263.3 

100.3 

0.2067 

0.00452 

998 

749.2 

0.0123 

-259.0 

102.3 

0.2067 

0.00385 

1087 

663.1 

0.0033 

-205.5 

103.1 

0.0400 

0.00426 

1119 

757.7 

0.0033 

-229.5 

103.9 

0.0400 

0.0043]| 

1139 

810.9 

0.0033 

-219.1 

102.5 

0.0400 

0. 00481 


.216 

.219 

.223 

.224 

0.215 

0.219 

0.229 

0.223 

0.241 

0.227 

0.247 

0.254 

0.259 


9.17 

12.6 

13.5 

13.7 

14.7 

14.5 

11.7 
11.3 

13.7 
13.0 

9.43 

11.5 

11.6 


87.8 1. 
85.2 1. 

83.7 1. 

84.1 1. 
87.4 1. 

84.9 1. 

80.8 1. 

81.9 1. 

80.9 1. 
82.8 1. 
76.7 1. 
77.0 1. 

73.2 1, 
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3.4.4 .1 Nominal Conditions-- Baseline Comparison. Hot— fire c* efficiency results 

for the four baseline elements are presented in Fig. 3-57. Values of p are 

c* 

shown as functions of the operating static pressure (not the computed throat stag- 
nation pressure) . The elements were throttled over a range of chamber pressure 
to ensure that a valid interpolation to the baseline 800 psia could be made. Data 
obtained at several values of center-post recess are presented for the baseline 
circular element and the baseline AR = 3 element. 

In general, the perforamnce dropped off slightly with increased chamber pressure 
(throttling up) . This is due to the fact that the gas velocity remains roughly 
constant while the liquid velocity increases, thus reducing V_/V . This is detri- 

L 

mental to both mixing and atomization. For the two elements fired with center-post 
recess, the performance was improved substantially as recess was increased. 

The variation of c* efficiency with aspect ratio, at zero recess and nominal condi- 
tions, is presented in Fig. 3-58. The performance first improved markedly with 
aspect ratio and then suddenly fell off again at the highest aspect ratio. 

The effect of center-post recess upon n c * is shown in Fig. 3-59 for the circular 
and the AR = 3 baseline elements at nominal conditions. For both elements, the 
performance was significantly improved with recess. However, the difference in 
performance between the elements was maintained, approximately. The performance 
of the recessed (R/y^ = 1.411) AR= 3 rectangular concentric tube element was 
notably high for a single-element injector (r| * = 95 percent). 

3.4.4. 2 Variations About AR = 3 . Characteristic velocity efficiency for the 
three variations about the aspect ratio 3.0 baseline element is presented in 
Fig. 3-60. Detailed comparison of the data with those obtained with the baseline 
AR = 3 element are presented in the next section in which the cold-flow and hot- 
fire results are correlated. One trend of note that is evident in these results 
is that the scaled down model of the baseline AR = 3 element (GCR-7) appears to 
be more sensitive to throttling than its larger counterpart. This is actually 

an indication that the smaller element is following the slopes predicted from 
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Figure 3-57. Characteristic Velocity Efficiency as a Function of 
Chamber Pressure (Baseline Elements) 
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Figure 3-58. 


Variation of c* Efficiency With Aspect Ratio (base- 
line elements, nominal conditions) 
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Figure 3-59. Characteristic Velocity Efficiency as a Function of 
Center-Post Recess (nominal conditions) 
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Figure 3-60. Variation About AR=3 Comparison 
Function of Chamber Pressure 
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cold flow more closely than the large element. The effect could be attributable 
to the interference of the chamber walls in the mixing and atomization process* 

The smaller element is "farther" from the walls and its flowfield would be affected 
to a lesser degree* 


3. S CORRELATION OF THE TEST RESULTS 


This section contains a comparison of the gas/ liquid cold-flow and hot-fire results. 

The approach to this comparison was discussed, briefly, in Section 3*1. Basically, 

the cold-flow results have been transformed into mixing limited Ol * . ) vapor- 

c mix 

ization limited (n * ) c* efficiency values and these values employed as the 


vap 

primary basis of comparison with hot-fire c* efficiency trends. 


The methods employed to compute the values of v\ * and v\ * . are summarized in 
r ' c vap c mix 

Appendix F. The results of these computations are presented in Fig. 3-61 and 3-62. 


The mixing limited c* efficiency, n * . , is presented in Fig. 3-61 as a function 

c mix 

of E |n for all of the injector elements of Phase III. It can be seen that one 

curve describes these data quite well. This curve, however, is only valid for 

mixture ratio 6.0. Other mixture ratios will yield separate curves. This points 

to the fact that E m and n * . are not uniquely related, 
m c mix 

The vaporization limited efficiency results presented in Fig. 3-62, wherein the 

value of n is shown as a function of the mass median liquid oxygen droplet 

c * 

diameter. The P droplet distribution function employed to compute these results 
was shown in Fig. 3-40. 

Comparisons between cold-flow and hot-fire are made at the nominal conditions. A 

list of these conditions is presented in Table 3-4. The results of the calculations 

of mixing limited and vaporixation limited efficiency are summarized in Table 3-1S. 

For each element, the values of E ffl , D wax > D L0X , V vap ’ ^ the P roduct of 

n * and n * (ri * „ ) are presented. 

c mix c vap c C.F. 
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PERCENT 




PERCENT (INVERSE LOG SCALE) 



0 100 200 300 1»00 


MASS MEDIAN LOX DROPLET DIAMETER, MICRONS 

Figure 3-62. Vaporization Limited c* Efficiency for LOX/GH- 
at Mixture Ratio 6.0:1 for P = 800 psia and 2 
Contraction Ratio 2.3:1 c 
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TABLE 3-15 RESULTS OF MIXING LIMITED AND VAPORIZATION LIMITED C* EFFICIENCY 


CALCULATIONS BASED UPON COLD-FLOW DATA (NOMINAL CONDITIONS) 


Element 
Code No. 

Element 

Type 

Relative 

Recess 

R/y 1 

E . 

Cl 

% 

V 

mix 

D 

wax 

microns 

°L0X 

microns 

V 

vap 

v 

C C.F. 

GCC-2 

Circle 

0 

47.0 

75.1 

395 

150 

97.3 

73.1 



■ 

49.0 

76.7 

342 

130 

98.0 

75.2 




52.8 

79.8 

310 

118 

98.4 

78.5 



1.344 

57.2 

83.3 

281 

107 

98.6 

82.1 

GCR-9 

AR = 1.5 

0 

57.2 

83.3 

399 

152 

97.2 

81.0 

GCR-8 

VO 

II 

CM 

0 

64.0 

88.0 

488 

186 

94.8 

83.4 




67.5 

90.2 

413 

157 

96.9 

87.4 




74.5 

93.4 

363 

138 

97.8 

91.3 



1.411 

81.0 

95.6 

314 

120 

98.3 

94.0 

GCR-10 

AR = 6 

0 

65.8 

89.2 

428 

163 

96.6 

86.2 

GCR-7 

Size 

0 

64.4 

88.3 

287 

109 

98.6 

87.1 

GCR-3 

Gas gap 

0 

49.5 

77.2 

321 

122 

98.3 

75.9 


















The conversion from wax droplet diameter to liquid oxygen droplet diameter was 
based upon a property correction for gas/ liquid atomization developed by Ingebo. 
This correction has been employed for gas/liquid concentric tube results in a 
FLOX/methane injector study at Rocketdyne (Ref. 7) . The value of the correction 
was computed with Eq. 3-17: 




where 


(2E) 

_p_LOX 



WAX 


0.381 


(3-17) 


o = surface tension 
y = viscosity 
p = density 


The physical properties of the various fluids employed in this program are pre- 
sented in Appendix G. 

Comparisons between the cold-flow and hot-fire c* efficiencies are presented in 

Fig. 3-63 through 3-66 for the baseline elements and two of the variations about 

AR = 3. In each figure curves for both the cold-flow c* efficiency (ri_* - 

^ C .F . 

n c *m- x x r, c* va p) and the mixing limited c* efficiency, nc* m i x » are presented, (i.e., 
curves for r lc* va p = 100 percent) . Examination of Fig . 3-63 through 3-66 will 
show that the curves representing mixing limited performance are much more 
representative of the hot-fire results than the curves that include a vaporization 
limited component. This result is quite reasonable since no account of secondary 
droplet breakup has been taken in the calculation of the liquid oxygen droplet 
diameters . 


In this thrust chamber (with a contraction ratio of 2.3. :1) the combustion gas 
can be expected to reach a velocity of 1400 ft/sec in the chamber and 5000 ft/sec 
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RELATIVE RECESS, R/y, 


Figure 3-64. Comparison of Cold-Flow and Hot-Fire c* Efficiencies 
for the CCTE (GCC-2) and the RCTE (AR=3, GCR-8) 
Presented as a Function of Center-Post Recess 
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Figure 3-65. Comparison of Cold-Flow and Hot-Fire c* 

Efficiencies Showing the Effect of Element 
Size 
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Figure 3-66. Comparison of Cold-Flow and Hot-Fire c* Efficiencies 
Showing the Effect of Gas Port Aspect Ratio, x^/y^ 
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at the throat. These velocities are more than adequate to produce substantial 
secondary droplet breakup. Rocketdyne is presently initiating a program (Ref. 8) 
to investigate the effect of secondary atomization with G/L circular concentric 
tube injector elements. 


No attempt was made to correct for the influence of secondary breakup in this 
report due to the fact that there were no experimental data available for G/L 
secondary breakup at the time of publication. Sufficient basic data have been 
presented in this document to allow corrections to be made at such time that these 
complementary data become available. However, the mixing limited c* results 
approximate the hot-fire data so well that, even with a more sophisticated droplet 


size correction, it appears that the values of n * 

c vap 

in most cases. 


must approach 100 percent 


The nominal conditions--baseline element comparison presented in Fig. 3-63 shows 
that three of the four elements follow the mixing limited c* trend almost exactly. 
The highest aspect ratio element does not fall on this curve. One possible explana- 
tion for this deviation is the following. As the aspect ratio is increased, rela- 
tive distortion of the flowfield increases (see Figure 3-29) . The aspect ratio 
6:1 element is the only element for which the liquid component significantly 
"breaks through'.' the gas component to the outer portions of the flowfield (see 
Fig. 3-29). This could mean that a significant portion of the liquid oxygen may 
have impinged upon the chamber wall forming a liquid film that would significantly 
interfere with both the mixing and vaporization mechanisms. This effect would 
only be associated with extremely small diameter chambers. In large chambers, 
the higher aspect ratio elements would be oriented with respect to the chamber 
wall in such a manner that the gas (fuel) component would be next to the wall, 
not the liquid component. 

A comparison of cold-flow and hot-fire c* efficiencies for the baseline CCTE and 
the AR = 3 RCTE is presented in Fig. 3-64 as a function of center-post recess. 

Here again, the mixing limited cold-flow curves are more representative of the 
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hot-fire performance. The level and increase of performance with recess are quite 
closely approximated in the case of the RCTE. The performance of the CCTE increased 
to a greater degree than was estimated from cold flow. 

Two of the variations about aspect ratio 3:1 are compared on the basis of c* effi- 
ciency in Fig. 3-65 and 3-66. In Fig. 3-65, the effect of element size clearly 
demonstrates the "mixing limited" nature of the hot-fire results. Reference to 
Fig. 3-32 and 3-44a will recall that the cold-flow results showed no change in 
level of mixing with size. However, a significant reduction of primary dropsize 
was suggested as element size was reduced. The hot-fire results show no significant 
change in performance with size. 

The effect of gas port aspect ratio is shown in Fig. 3-66. The trends suggested 
by the cold-flow results are drawn only to represent the change from one aspect 
ratio to the other, and are not representative of the true curves that should 
connect these two points (see Fig. 3-35). The change in level of performance, 
and indeed the actual value of performance, are represented by the mixing limited 
cold-flow results. 

In summary, it may be concluded that the correlation between the mixing limited 
cold-flow c* efficiencies and the hot-fire performance values is excellent. It 
may be assumed, therefore, that the mixing trends suggested by the cold-flow data 
are truly representative of the processes encountered in hot fire and that the cold- 
flow design guidelines may be employed in the formulation of injector designs for 
hot fire. 

Due to the mixing limited nature of the hot-fire results, the trends in concentric 
tube atomization characteristics have neither been confirmed nor denied. Other 
studies, however, have clearly confirmed the validity of the cold-flow atomization 
representation (Ref. 7). 
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4.0 LIQUID/LIQUID STUDIES 


Although Phases III and IV of the subject program were primarily dedicated to the 
investigation of gas/liquid injector elements, two significant liquid/liquid in- 
jector studies were undertaken. These studies were conducted during Phase IV, 
and included a mixing and atomization study of unlike triplet elements with cir- 
cular orifices and a mixing study of unlike-doublet elements with rectangular 
orifices having rounded entrances. 

The triplet studies were conducted to gain a better appreciation of the mixing 
and atomization levels achievable with triplet elements. This was to serve as a 
basis for comparison of the performance of triplet injector elements with the 
performance of unlike-doublet elements. 

The unlike-doublet studies were conducted to obtain mixing data for rectangular 
unlike-doublet elements having orifices with rounded inlets. These results com- 
plement the rectangular unlike-doublet mixing results of Phase I which were ob- 
tained with sharp entranced orifices. 

No attempt has been made in this document to give the background for the liquid/ 
liquid studies nor the interrelationship and correlation of these data with other 
liquid/liquid results obtained during Phase I. A detailed discussion of all of 
the liquid/liquid results is presented in the Comprehensive Program Summary Re- 
port (Ref. 3). A full appreciation for these data and their implication for in- 
jector design can only be obtained by reading Ref. 3. The results obtained dur- 
ing Phase IV are presented in this report for completeness of documentation, only. 
It was decided that a comprehensive discussion of the liquid/liquid results in 
both the Program Summary and this document would be unwarranted. 

The mixing results for both the triplets and the unlike doublets are presented in 
terms of the "centerline momentum ratio," <j>. (See Appendix H and Ref. 3 for back- 
ground.) The physical basis for the parameter is identical for the different ele- 
ment types. However, the computation of this parameter is dependent upon the 
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element type. The equations for the definition of the centerline momentum ratio 
are presented below for triplet and uni ike-doublet elements: 


• triplet 



»2 V 2 V 2 


(4-1) 


• rectangular unlike doublet 


(4-2) 


In the case of the triplet elements, the outer orifices are denoted by the sub- 
script 1 and the central orifice by the subscript 2. For the unlike doublets, 
the larger orifice (area) is denoted by the subscript o (oxidizer) and the smaller 
orifice by the subscript f (fuel). For this particular study, the widths of the 
o and the f orifices were the same. 

4.1 TEST HARDWARE 

The design specifications for circular triplet elements and the rectangular 
unlike-doublet elements are presented in Table 4-1. The nomenclature employed 
here are defined in Fig. 4-1. 


<J> = 


. p f V f b f 


p V 2 b 
0 o o 
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TABLE 4-1. 


NONCIRCULAR ORIFICE PROGRAM- -PHASE IV LIQUID/LIQUID STUDIES HARDWARE SUMMARY 


Element , 
Code °1 

No. (in.) 

d 2 

(in.) 

V d 2 

b 

0 

(in.) 

b f 

(in.) 

V b o 

w 

(in.) 

b o /w 

AR 

0 

Al 2 
(in. ) 

A 2 
4 2 
(in. ) 

A 

o 2 

(in.) 

Af 2 
(in. ) 

th 

CD i 

& 

4 


UT-1 

0.0300 

0.0300 

1.000 

mm—w. 

_ _ __ 

_ _ 

_ _ „ _ 

... 

0.001414 

0.000707 





13.5 

15.0 

m 

UT-2 

0.0300 

0.0430 

0.698 

— 

— 

— 

— 

— 

0.001414 

0.00145 

— 

— 

13.5 

10,5 

-- 

UT-3 

0.0300 

0.0575 

0.522 

— 

— 




— 

0.001414 

0.00260 

— 

— 

13.5 

7.83 

— 

UT-4 

0.0200 

0.0200 

1.000 

— 

— 

— 

---- 

— 

0.000628 

0.000314 

— 

— 

13.5 

15.0 

— 

ITT-S* 

0.0670 

0.0670 

1.000 

— 

— 

— 

— 

— 

0.00706 

0.00353 

— 

---- 

50.0 

50.0 

— 


UD-1 

— 

— 

— 

0.0490 

0.0490 

1.000 

0,0490 

1.000 




0.00240 

0.00240 

-- 

— 

10 

UD-2 

— 

— 

— 

0.0930 

0.0930 

1.000 

0.0258 

3.600 

— 

— 

0.00240 

0.00240 


-- 



UD-3 

— 

— 

— 

0.1200 

0.1200 

1.000 

0.0200 

6.000 

— 

— 

0.00240 

0.00240 

-- 

-- 



UD-4 

— 

— 

— 

0.0327 

0.0204 

0.625 

0.1175 

0.278 

— 

— 

0.00384 

0.00240 

-- 

— 



UD-5 

— 

— 

— 

0.0620 

0.0388 

0.625 

0.0620 

1.000 

— 

— 

0.00384 

0,00240 

-- 

— 



UD-6 

— 

— 

— 

0.0736 

0.0736 

0.625 

0.0327 

3.600 


— 

0.00384 

0.00240 

-- 




UD-7 

— 

— 

— 

0.1518 

0,0949 

0.625 

0.02S3 

6.000 

— 

— 

0.00384 

0.00240 


-- 



UD-8 

— 



— 

0.0783 

0.0306 

0.391 

0.0783 

1.000 

— 

.... 

0.00614 

0.00240 

-- 

— 



UD-9 

— 

— 

— 

0.1487 

0.0581 

0.391 

0.0413 

3.600 

— 



0.00614 

0.00240 

— 




UD-10 

— 

— 

— 

0.1920 

0.0751 

0.391 

0.0320 

6.000 

— 



0.00614 

0.00240 

-- 


10 



> Triplei 


10 1 


> Unlike 


10 


♦Results for 


this element were obtained under NASA Contract 


No. NAS7-726 and 


are included here for completeness. 


Doublets 



• TRIPLET ELEMENTS 



• UNLIKE DOUBLETS 


LARGER 


ORIFICE 


SMALLER 

ORIFICE 


KH h b .H 




Figure 4-1. Triplet and Unlike-Doublet Nomenclature 
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The model hardware assemblies for the two element types were quite similar. A 
schematic cross section of a typical triplet model is presented in Fig. 4-2 and 
that for a typical unlike-doublet model in Fig. 4-3. 

The triplet model configuration (Pig* 4-2) is typical for the two elements having 
a diameter ratio unity, d^/d^ = 1.0. The elements are labeled UT-1 and UT-4. The 
diameter ratios 1.433 (UT-2) and 1.917 (UT-3) were obtained by drilling out the 
central orifice of element UT-1 (see Table 4-1) . Element UT-5 was not fabricated 
or tested under this contract. This element was tested under NASA contract 
NAS7-726, and is documented in Ref. 9. Atomization results from UT-5 were incor- 
porated in the discussion of triplet atomization in this report for completeness. 
All triplet elements incorporated circular orifices with rounded inlets. The outer 
orifices were canted toward the central orifice with an included angle of 60 
degrees. 

The rectangular unlike doublets (Fig. 4-3) were configured with a face orifice 
spacing of 0.26 inch, rounded inlets, orifice length-to-hydraulic diameter ratio 
10:1, and an included impingement angle of 60 degrees. 

The model element injector body design was similar for the triplet and doublet 
elements. The only difference being the number of orifice inlet fittings — two 
for the doublets and three for the triplets. A typical injector body is shown 
in Fig. 4-4 (unlike doublet shown) . 

Face patterns for three uni ike- doublet models (UD-1, UD-2, UD-3 - see Table 4-1) 
are shown in Fig. 4-5. The orifices for the rectangular uni ike-doublet models 
were fabricated by the electrical discharge machining (EDM) method. 

4.2 EXPERIMENTAL RESULTS 

The liquid/liquid experimental results are presented in two parts; triplet re- 
sults and unlike-doublet results. The triplet results discussion has a mixing 
and an atomization section, while the discussion of unlike-doublet results has 
only a mixing section. 
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ORIFICE NO. 2 



Figure 4-2. Schematic of Triplet Element Geometry (Typical for Models With 
d 2 /d. = 1, see Table 4-1) 




0 . 26 " 


Figure 4-3. Rectangular Uni ike -Doublet Model Format Typical for 
all 10 Models 
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5AA34-5/21/73-S1A 


Figure 4-4. Model Injector Element Body 
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Figure 4-5. Typical Rectangular Unlike Doublet Face Patterns 


4.2.1 Triplet Results (Circular Orifices) 


4. 2. 1.1 Triplet Mixing Results . The triplet mixture ratio uniformity results 
are presented in Table 4-2. Reference may be made to Table 4-1 for model con- 
figuration descriptions. The results are interpreted in terms of the mixture 
ratio uniformity parameter, (see Appendix F), and the centerline momentum 
ratio, <{> (see Eq. 4-2 and Appendix H) . 

The effects of orifice diameter, diameter ratio, and centerline momentum upon 
the mixing of the triplet elements tested are shown in Fig. 4-6. It is quite 
evident that the mixing level achieves a maximum value for values of $ near 
unity. Actually, for the elements having d^A^ s 1.00, the optimum E m was 
achieved for <f> = 1.00, while for the elements having d^/d^ < 1.0, the optimum 
values of E m were obtained at a value of <J> = 0.714. This suggests that, for 
triplets, the specification of both 4> and required to ensure optimum 

mixing. That is to say: 

E m = d l /d 2^ (4 “ 3) 

opt 

The effect of orifice diameter is shown in Fig. 4-6a for dj/d 2 = 1.00. The re- 
sults were obtained with elements having 0.020- and 0.030-inch diameter orifices. 
It can be seen that at least for this magnitude of size variation, E m is rela- 
tively insensitive to the orifice diameter, especially at the optimum point, 

= 1 . 0 . 

The effect of diameter ratio is shown in Fig. 4-6b. Two effects are evident: 

(1) the value of 4> for optimum mixing moves with diameter ratio, and (2) the 
level of mixing drops with decreasing diameter ratio, d^/d^. This latter in- 
fluence is highlighted in the cross plot of the data presented in Fig. 4-6c. 

(The dashed curve in Fig. 4-6 represents a correction of the results to equal 
total flowrate- -correction obtained from Fig. 4-7.) Here, the maximum values 
of E m are plotted with respect to diameter ratio. (Note: Values of d^/d^ < 1.0 

represent triplet elements with the central orifice larger than the outer 
orifices . ) 
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TABLE 4-2. TRIPLET MIXING RESULTS 


TEST 

NO. 

ELEMENT 

CODE 

NO. 

"t 

(lbm/sec) 

*2 

(lbm/sec) 

"i 

(lbm/sec) 

Vi 

V 2 

V 1 

4> 

Em,% 

1 

UT-1 

0.0777 

0.0348 

0.0429 

0.811 

77.7 

70.0 

1.109 


2 



0.0273 

0.0504 

0.542 

70.0 

82.2 

2.483 


3 



0.0427 

0.0350 

1.220 

95.4 

57.1 

0.490 

83.3 

4 



0.0502 

0.0275 

1.825 

112.0 

44.9 

0.219 

69.6 

5 

UT-1 

0.0777 

0.0206 

0.0571 

0.361 

46.0 

93.2 

5.596 

45.2 

6 

UT-4 

0.0345 

0.0154 

0.0191 

0.811 

77.7 

70.0 

1.109 

91.9 

7 



0.0121 

0.0224 

0.542 

61.0 

82.2 

2.483 

71.4 

8 

UT-4 

0.0345 

0.0190 

0.0155 

1.220 

95.3 

57.1 

0.490 

79.4 

9 

UT-2 

0.108 

0.0652 

0.0428 

1.522 

71.0 

69.9 

0.917 

86.5 

10 



0.0751 

0.0329 

2.283 

81.8 

53.7 

0.408 

81.2 

11 



0.0544 

0.0536 

1.015 

59.2 

87.5 

2.062 

54.6 

12 



0.0835 

0.0245 

3.416 

91.0 

39.9 

0.182 

58.8 

13 



0.0435 

0.0645 

0.675 

47.4 

105.0 

4.660 

28.0 

14 



0.0681 

0.0399 

1.709 

74.2 

65.1 

0.727 

90.4 

15 



0.0652 

0.0428 

1.522 

71.0 

69.9 

0.917 

88.4 

16 


0.144 

0.0908 

0.0532 

1.709 

98.9 

86.8 

0.727 

87.3 

18 

UT-2 

0.0811 

0.0512 

0.0299 

1.709 

SS.7 

48.9 

0.727 

90.7 

19 

UT-3 

0.165 

0.1251 

0.0399 

3.133 

76.0 

65.2 

0.523 

80.0 

20 



0.1187 

0.0463 

2.562 

72.1 

75.6 

0.782 

82.3 

21 



0.1213 

0.0437 

2.773 

73.6 

71.4 

0.668 

82.2 

22 



0.1136 

0.0514 

2.208 

69.0 

83.9 

1.053 

81.0 

23 



0.1301 

0.0349 

3.725 

79.0 

57.0 

0.360 

71.2 

24 

UT-3 

0.165 

0.1074 

0.0576 

1.863 

65.2 

94.1 

1.479 

78.6 



























0 0.05 0.10 0.15 0.20 

TOTAL ELEMENT FLOWRATE, LBM/SEC 


Figure 4-7. Effect of Total Flowrate Level Upon Unlike Triplet 
Mixing 
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The effect of total element flowrate upon triplet mixture ratio uniformity is 
shown in Fig. 4-7. A 4 to 5 percentage point reduction in E was encountered 
for a factor of two increase in total flow. 

4,2, 1.2 Triplet Atomization Results . The results of the triplet atomization ex- 
periments are summarized in Table 4-3 for the outer orifices and Table 4-4 for 
the central orifices. All of the atomization tests were conducted at the value 
of centerline momentum ratio, <J>, which was found to yield optimum mixing for each 
individual element . 

Results in Tables 4-3 and 4-4 are presented in terms of the mass median droplet 
diameter, D, and the relative droplet diameter, D/Dj (droplet diameter to orifice 
diameter ratio). 

A detailed analysis of these results suggested that only the dropsize data ob- 
tained for wax jet Reynolds numbers less than 10 4 (laminar jets) could be corre- 
lated. The results for Reynolds numbers above this value indicated a great deal 
of scatter due to the fact that the jets were in transition from laminar to tur- 
bulent flow (Ref. 3 and 4) . 

The best correlation of the laminar jet atomization data was found to be of the 
form: 

D/ D j = f (W e , L/D, R e ) (4-4) 

where 

W = Weber number = pV^D./cr 
e J )' 

R = Reynolds number = pV.D./y 

J J 

L/D = orifice length to diameter ratio 
Vj = wax jet velocity 

p = wax density 

p = wax viscosity 

o = wax surface tension 
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TABLE 4-3. TRIPLET ATOMIZATION RESULTS, WAX IN OUTER ORIFICES, NO. 1 


TEST 

NO. 

ELEMENT 

CODE 

NO. 

w T 

(Ibm/sec) 

*2 

(lbm/sec) 

(lbm/sec) 

V»i 

WAX 

REYNOLDS 

NO. 

WAX 

VELOCITY 

(ft/sec) 

D, u 

D/D j 

1 

UT-4 

0.01080 

0.00465 

0.00615 

0.756 

1175 

30 

640 

1.260 

2 


0.01794 

0.00774 

0.0102 

0.759 

1959 

50 

341 


3 


0.0360 

0.0155 

0.0205 

0.756 

3917 


205 


4 


0.0541 

0.0233 

0.0308 

0.756 

5876 

150 

157 

0.309 

5 

UT-4 

0.0720 

0.0310 

0.0410 

0.756 

7834 


123 

0.242 

6 

UT-1 

0.0243 

0.0105 

0.0138 

0.761 

1763 

30 

620 

0.814 

7 


0.0405 

0.0174 

0.0231 

0.753 

2938 

50 

428 

0.562 

8 


0.0810 

0.0349 

0.0461 

0.757 

5876 

100 

240 

0.315 

9 


0.1215 

0.0523 

0.0692 

0.756 

8813 

150 

136 

0.178 

10 

UT-1 

0.1621 

0.0698 

0.0923 

0.756 


187 

148 

0.194 

29 

UT-2 

0.0610 

0.0379 

0.0231 

1.641 

2938 

50 

515 

0.676 

30 


0.1218 

0.0757 

0.0461 

1.642 

5876 

100 

291 

0.382 

31 


0.1829 

0.1137 

0.0692 

1.643 

8813 

150 

197 

0.2S9 

32 


0.2438 

0.1515 

0.0923 

1.641 

11752 

200 

175 

0.230 

33 

UT-2 

0.3657 

0.2273 

0.1384 

1.642 

17628 

300 

175 

0.230 

34 

UT-3 

0.0819 

0.0588 

0.0231 

2.545 

2938 

50 

490 

0.643 

35 


0.1638 

0.1177 

0.0461 

2.553 ! 

5876 

100 

200 

0.263 

36 

UT-3 

0.2457 

0.1765 

0.0692 

2.551 1 

8813 

150 

200 

0.263 




DATA FROM NAS 7- 

•726 (Ref. 9) 




185 

UT-5 

0.1081 

0.0457 


0.732 

3504 

27 

583 

0.343 

186 


0.1636 

0.0700 

0.0936 

0.748 

5250 

40 

470 

0.276 

187 


0.2186 

0.0943 

0.1243 

0.759 

6969 

53 

412 

0.242 

188 


0.3269 

0.1408 

0.1861 

0.757 

10434 

80 

307 

0.180 

189 


0.4292 

0.1853 

0.2439 

0.760 

13689 

104 

266 

0.156 

190 


0.5603 

0.2393 


1 . 

18007 

137 

190 

0.112 

191 

UT-5 

0.6123 

0.2583 

0.3540 


19857 

151 

195 

0.115 
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TABLE 4-4. TRIPLET ATOMIZATION RESULTS, WAX IN CENTRAL ORIFICE, NO. 2 


TEST 

NO. 

ELEMENT 

CODE 

NO. 

“t 

(lbm/sec) 

*2 

(lbm/sec) 

(lbm/sec) 

v«i 

WAX 

REYNOLDS 

NO. 

WAX 

VELOCITY 

(ft/sec) 

D, 

'V 

D/D j 

12 

UT-1 

0.01852 

0.00692 

0.0116 

0.597 

1763 

30 

820 

wm 

13 


0.0309 

0.0115 

0.0194 

0.593 

2938 

50 

405 

mmm 

14 


0.0618 

0.0231 

0.0387 

0.597 

5876 


208 

9 UK 

15 


0.0927 

0.0346 

0.0581 

0.596 

8813 


155 

0.203 

16 


0.1235 

0.0461 

0.0774 

0.596 

11751 


108 

0.142 

17 

UT-1 

0.1853 

0.0692 

0.1161 

0.596 

17627 


96 

0.126 

18 

UT-4 

0.00825 

0.00308 

0.00517 

0.596 

1175 

30 

840 

1.654 

19 


0.01373 

0.00513 

0.00860 

0.597 

1959 

50 

355 

0.699 

20 


0.0275 

0.0103 

0.0172 

0.599 

3917 

100 

172 

0.339 

21 


0.0412 

0.0154 

0.0258 

0.597 

5876 

150 

140 

0.276 

22 


0.0549 

0.0205 

0.0344 

0.596 

7834 

200 

114 

0.224 

23 


0.0824 

0.0308 

0.0516 

0.597 

11751 

300 

112 

0.221 

24 


0.0433 

0.0241 

0.0192 

1.255 

4212 

50 

611 

0.559 

25 


0.0864 

0.0481 

0.0383 

1.256 

8422 

100 

307 

0.281 

26 


0.1297 

0.0722 

0.0575 

1.256 

12633 

150 

190 

0.174 

27 


0.1728 

0.0962 

0.0766 

1.256 

16843 

200 

154 

0.141 

28 

UT-2 

0.2593 

0.1443 

0.1150 

1.255 

25265 

300 

149 

0.136 

39 

UT-3 

0.0643 

0.0424 

0.0219 

1.936 

5631 

50 

800 

0.548 

40 


0.1286 

0.0848 

0.0438 

1.936 

11262 

100 

372 

0.255 

41 


0.1929 

0.1271 

0.0658 

1.932 

16892 

ISO 

209 

0.143 

42 


0.2572 

0.1695 

0.0877 

1.933 

22523 

200 

220 

0.151 

43 

UT-3 

0.3859 

0.2543 

0.1316 

1.932 

33777 

300 

173 

0.119 




DATA FROM NAS 7 

-726 (Ref. 9) 




192 

UT-5 

0.0843 

0.0315 

0.0528 

0.597 

3530 

27 

606 

0.356 

193 


0.1258 

0.0465 

0.0793 

0.586 

5223 

40 

484 

0.284 

194 


0.1650 

0.0629 

0.1021 

0.616 

7060 

54 

425 

0.250 

195 


0.2477 

0.0935 

0.1542 

0.606 

10498 

80 

333 

0.196 

196 


0.3294 

0.1201 

0.2093 

0.574 

13477 

103 

297 

0.175 

197 

UT-5 

0.4155 

0.1572 

0.2583 

0.609 

17636 

134 

238 

0.140 














(It must be noted that the values of p, y, and a were not varied and thus the 
general specification of Weber number and Reynolds number as independent param- 
eters is speculative.) 

Equation 4-4 states that dropsize for triplet elements in laminar flow is a func- 
tion of Weber number and orifice L/D. The Weber number describes the dynamic as- 
pects of the atomization process while the L/D and Reynolds number describe the 
degree of jet velocity profile development for laminar flow in orifices with 
rounded entrances. 

Atomization results are presented in Fig. 4-8 for the outer orifices and in Fig. 
4-9 for the central orifices. The dropsizes are presented in terms of Weber num- 
ber and L/D only, while the effect of Reynolds number has been neglected. The 
following argument is presented to show that the Reynolds number effect is truly 
second order, at least over the range of variables tested. 

Due to the fact that the physical properties of the fluids (p, y, a) were not 
varied, the Reynolds number, Weber number, and jet diameter are related, uniquely: 

R = constant x J W x D . (4-5) 

0 6 J 

Thus, for a constant jet diameter, the independent effect of Reynolds number and 
Weber number cannot be determined. However, the results shown in Fig. 4-8 and 
4-9 suggest that as long as the orifice L/D is equivalent, the data obtained with 
different jet diameters are well represented by a single curve. This strongly 
suggests that the independent effect of Reynolds number is, indeed, second order. 

It must be emphasized that these representations do not imply the independent in- 
fluence of fluid density, viscosity, or surface tension. Also, no account is 
made of the dynamic condition or physical properties of the jets opposing the wax 
jets. The data were obtained with hot water as a propellant simulant injected in 
those orifices opposing the wax (see Appendix D) . 
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An interesting cross correlation of the data is presented in Fig. 4-10, wherein 
[)/Dj is plotted as a function of the product of Weber number and orifice L/D for 
both central and outer jets. A single curve is representative of the data for 
both No. 1 and No. 2 orifices in the range 10 5 < (W ) (L/D) < 3 x 10 6 . For values 

c ® 

of (W g ) (L/D) below 10°, the curves diverge. An analytical expression of the 
curve in the range 10^ to 3 x 10^ is: 


D = 2.15 


x 10 5 V.' 0 * 77 (L/D)“ 0,39 D. 0 * 62 
3 3 


(4-6) 


where 


V. = ft/sec 
Dj = inch 
D = microns 

A comparison of the atomization characteristics (laminar wax jets) of an unlike- 
triplet element, an unlike-doublet element and a like-doublet element is pre- 
sented in Fig. 4-11. The comparison is made for the fuel orifice of each ele- 
ment. The orifice diameter is 0.067 inch and the L/D = 50. Dropsize values for 
all elements were obtained from correlations presented by Zajac (Ref. 9) for lam- 
inar jets. The diameter ratios for the triplet, the like doublet, and the unlike 
doublet are unity. 

The like and unlike doublets produce approximately the same dropsize while the 
triplet produces somewhat larger sizes. 

4.2.2 Unlike-Douhlet Results (Rectangular Orifices) 


The objective of the unlike-doublet experiments conducted during Phase IV was to 
obtain critical design data for rectangular doublets with rounded entrances to 
complement the data formerly obtained during Phase I for doublets with sharp en- 
trances. Primary experimental effort was centered about the parametric investi- 
gation of the effect of element area ratio and orifice aspect ratio upon mixing. 
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Figure 4-10. Correlation of Relative Droplet Diameter With Weber No. and Orifice L/D for 
Unlike Triplets (Both Central and Outer Orifice Results Shown) 















MICRONS 



WAX INJECTION VELOCITY, FT/SEC 


Figure 4-11. Comparison of Wax Droplet Diameters for a Triplet, 
an Unlike Doublet, and a Like Doublet (Laminar 
Wax) 
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The geometrical and operational parameters that were investigated are presented in 
Fig. 4-12. 


Centerline momentum ratio, <£ (see Appendix H) , has been shown to be the most criti- 
cal correlating parameter for unlike-doublet mixing. The mixing level, E^, of 
single elements, optimizes normally at a value of <j> % 1.0. As a matter of fact, 
the expression <J> = 1 has, for years, served as a design criterion for unlike- 
doublet injectors. When a propellant combination is specified (i.e., P Q /Pf) at a 
given mixture ratio and the condition w q /w£ = 1 is imposed, the value of b Q /b£ 
is defined for the case <f> = 1. Thus, the only independent variable remaining is 
the aspect ratio of one of the orifices (e.g., b 0 /w Q = AR q ) . The aspect ratio of 
the other orifice is functionally related to AR q and is not an independent 
parameter. 

The restriction of equal facing widths, w q /w£ = 1 (Fig. 4-12) was imposed to 
limit the scope of the investigation. A logical extension of the work reported in 
this document would be the investigation of elements for which w Q /w£ f 1. 

The test matrix that was selected for the unlike-doublet experiments is shown in 
qualitative form in Fig, 4-13. Three values of area ratio A^/A^ were selected 
(1.00, 0.625, and 0.391). (Note: due to the restriction w q = w f , A Q /A f = b o /b f .) 

For each value of area ratio, b Q /b£, at least three values of oxidizer orifice as- 
pect ratio b 0 /w Q were selected (four values were selected for b Q /b£ = 0.625). The 
element design specifications are summarized in Table 4-1. 

4.2.2. 1 Unlike-Doublet Mixing Results . The rectangular unlike-doublet, cold-flow 
mixing data are presented in Table 4-5. The results are interpreted in terms of 
the mixture ratio uniformity parameter, E m (see Appendix F) , as a function of the 
centerline momentum ratio parameter, <p (see Appendix H) . 

The basic mixing characteristics of the 10 unlike-doublet elements are presented 
in Fig. 4-14 through 4-16. Each figure is reserved for results obtained at one 
value of area ratio, t>£/b Q . Within each figure, the data for the various oxidizer 
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Element Nomenclature 


b f Fuel 


Basic Definitions and Ground Rules 

w = w £ = w 
o f 

AR = b /w 
o o 

AR f = b f /w 

Area Ratio = A.p/A Q = b^/b D 

p /p £ = const = 1 .44 
o f 

6 = impingement angle = 60 degrees 
Correlating Parameter 

Centerline Momentum Ratio 


4 > = 


°f b f 



b 

o 




where 

MR = w /w, 
o T 


Figure 4-12. Definition of Nomenclature for Unlike-Doublet 
Cold-Flow Testing 
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TABLE 4-5. RECTANGULAR UNLIKE-DOUBLET MIXING RESULTS, ROUNDED ENTRANCES 



ELEMENT 

CODE 

NO. 


(lbm/sec) (lbm/sec) 



0.099 

0.118 

0.131 

0.083 

0.091 

0.099 

0.106 

0.114 

0.099 

0.101 

0.083 

0.091 

0.099 

0.106 

0.114 

0.111 

0.109 

0.128 

0.105 

0.120 

0.078 

0.094 

0.110 

0.130 

0.142 

0.094 

0.104 

0.109 

0.116 

0.123 


"f 

(lbm/sec) 

V K f 

V 

0 

(ft/sec) 

V f 

(ft/sec) 


0.086 

1.09 

61.8 

82.7 

1.25 

0.112 

0.607 

44.7 

107.0 

3.99 

0.097 

0.856 

54.6 

93.3 

2.01 

0.081 

1.22 

65.1 

77.9 

0.984 

0.062 

1.90 

77.6 

59.6 

0.406 

0.049 

2.67 

86.2 

47.1 

0.206 

0.097 

0.856 

54.6 

93.3 

2.00 

0.089 

1.02 

59.9 

85.6 

1.41 

0.081 

1.22 

65.1 

77.9 

1.00 

0.074 

1.43 

69.7 

71.2 

0.707 

0.066 

1.73 

75.0 

63.5 

0.500 

0.081 

1.22 

65.1 

77.9 

1.00 

0.079 

1.28 

66.4 

76.0 

0-.900 

0.097 

0.856 

54.6 

93.3 

2.00 

0.089 

1.02 

59.9 

85.6 

1.41 

0.081 

1.22 

65.1 

77.9 

1.00 

0.074 

1.43 

69.7 

71.2 

0.707 

0.066 

1.73 

75.0 

63.5 

0.500 

0.069 

1.61 

45.6 

66.3 

0.900 

0.071 

1.54 

44.8 

68.3 

1.00 

0.052 

2.46 

52.6 

50.0 

0.380 

0.075 

1.40 

43.2 

72.1 

1.18 

0.060 

2.00 

49.3 

57.7 

0.580 

0.102 

0.765 

.32.1 

98.1 

4.14 

0.086 

1.09 

38.7 

82.7 

1.93 

0.070 

1.57 

45.2 

67.3 

0.934 

0.050 

2.60 

53.5 

48.1 

0.341 

0.038 

3.74 

57.6 

36.5 

0.164 

0.086 

1.09 

38.7 

82.7 

1.97 

0.079 

1.32 

42.8 

76.0 

1.36 

0.071 

1.54 

44.8 

68.3 

0.998 

0.064 

1.81 

47.7 

61.5 

0.716 

0.057 

2.16 

50.6 

54.8 

0.505 

1 . 



78.5 

70.6 
61.2 
75.0 

85.7 

83.4 

75.0 

61.1 

86.8 

84.4 

72.4 
84.2 
88.8 
80.0 
68.8 

79.0 

79.6 

77.1 

75.8 

74.0 

79.8 

77.1 

78.9 

68.2 

69.7 

80.4 
85.0 

85.8 
80.2 

67.8 
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TABLE 4-5. (Continued) 


TEST 

NO. 

ELEMENT 

CODE 

NO. 

W T 

(lbm/sec) 

• 

W 

o 

(lbm/sec) 

(lbm/sec) 

W o /Wf 

V 

o 

(ft/sec) 

V f 

(ft/sec) 


Em, % 

22 

UD-6 

0.180 

0.117 

0.063 

1.86 

48.1 

60.6 

0.670 

68.4 

21 



0.102 

0.078 

1.31 

41.9 

75.0 

1.37 

82.6 

20 



0.094 

0.086 

1.09 

38.7 

82.7 

1.97 

70.6 

19 



0.130 

0.050 

2.60 

53.5 

48.1 

0.348 

48.6 

18 



0.110 

0.070 

1.57 

45.2 

67.3 

0.950 

86.6 

17 



0.110 

0.070 

1.57 

45.2 

67.3 

0.950 

87.9 

16 



0.110 

0.070 

1.57 

45.2 

67.3 

0.950 

79.1 

15 



0.078 

0.102 

0.765 

32.1 

98.1 

4.02 

42.8 

14 



0.094 

0.086 

1.09 

38.7 

82.7 

1.97 

61.2 

13 



0.110 

0.070 

1.57 

45.2 

67.3 

0.953 

86.2 

12 



0.130 

0.050 

2.60 

53.5 

48.1 

0.348 

53.0 

11 

UD-6 

0.180 

0.142 

0.038 

3.74 

58.4 

36.5 

0.169 

36.0 

38 

UD-7 

0.180 

0.109 

0.071 

1.54 

44.8 

68.3 

1.00 

82.9 

37 



0.111 

0,069 

1.61 

45.6 

66.3 

0.900 

86.1 

36 



0.114 

0.066 

1.73 

46.9 

63.5 

0.800 

85.2 

35 



0.126 

0.054 

2.33 

51.8 

51.9 

0.432 

63.9 

34 



0.130 

0.050 

2.60 

53.5 

48.1 

0.348 

55.9 

33 



0.134 

0.046 

2.91 

55.1 

44.2 

0.277 

50.0 

32 



0.094 

0.086 

1.09 

38.7 

82.7 

2.00 

57.4 

31 



0.101 

0.079 

1.28 

41.5 

76.0 

1.41 

68.0 

30 



0.109 

0.071 

1.54 

44.8 

68.3 

1.00 

78.8 

29 



0.116 

0.064 

1.82 

47.7 

61.5 

0.707 

85.3 

28 

UD-7 

0.180 

0.123 

0.057 

2.16 

50.6 

54.8 

0.500 

76.9 

70 

UD-8 

0.180 

0.123 

0.057 

2.16 

31.6 

54.8 

0.800 

86.9 

43 



0.104 

0.076 

1.37 

26.7 

73.1 

2.00 

69.3 

42 



0.112 

0.068 

1.65 

28.8 

65.4 

1.41 

75.2 

41 



0.119 

0.061 

1.95 

30.6 

58.7 

1.00 

82.8 

40 



0.126 

0.054 

2.33 

32.4 

51.9 

0.707 

84.7 

39 

UD-8 

0.180 

0.132 

0.048 

2.75 

33.9 

46.2 

0.500 

72.2 
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TABLE 4-5. (Concluded) 


TEST 

NO. 

ELEMENT 

CODE 

NO. 

w T 

(lbm/sec) 

i 

(lbm/sec) 

w f 

(lbm/sec) 

*<A 

V 

0 

(ft/sec) 

V f 

(ft/sec) 

4> 

Em, % 

74 

UD-9 

0.180 

0.119 


1.95 

30.6 

58.7 

1.00 

81.9 

71 



0.112 


1.65 

28.8 

65.4 

1.36 

83.4 

48 



0.104 

lill" 

1.37 

26.7 

73.1 

2.00 

81.1 

47 



0.112 

0.068 

1.65 

28.8 

65.4 

1.36 

74.2 

46 



0.119 

0.061 

1.95 

30.6 

58.7 

1.00 

83.7 

45 



0.126 

0.054 

2.33 

32.4 

51.9 

0.680 

85.8 

44 

UD-9 

0.180 

0.132 


2.75 

33.9 

46.2 

0.500 

74.4 

53 

UD-10 

0.180 

0.104 


1.37 

26.7 

73.1 

2.00 

72.2 

52 



0,112 


1.65 

28.8 

65.4 

1.41 

82.1 

51 



0.119 

1 mm '^heBe 

1.75 

30.6 

65.4 

1.00 

80.5 

50 



0.126 

0.054 

2.33 

32.4 

51.9 

0.707 

69.3 

49 

UD-10 

0.180 

0.132 

0.048 

2.75 

33.9 

46.2 

0.500 

54.1 
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Figure 4-1S. Rectangular Unlike-Doublet Mixing, Ejj, as a Function of 0 
and AR for b~/b = 0.625 (Rounded Entrances) 

















orifice aspect ratios are presented in separate plots. In all cases, it may be 
noted that the mixing level for a given element optimized at a value of <f> near 
unity. The value of the maximum level of mixing varies with both b f /b o and b Q /w, 
however. 

A cross plot of the data is presented in Fig. 4-17, wherein the maximum E m is 
shown as a function of AR q along lines of constant b^/b^ It is clear from the 
results that the optimum value of AR q is functionally related to the value of 
b^/b^. A correlation of all of the data is presented in Fig. 4-18, which sug- 
gests the analytical form of this relationship. Here, the maximum values of E 
are plotted as a function of AR Q /(b ;f /b 0 )' 5 . This presentation allows a single 
curve to be drawn through the data. Mixing results previously obtained with 
rectangular unlike doublets with sharp-edged entrances (see Ref. 1 and 3) are 
also presented in Fig. 4-18. For both rounded and sharp-edged entrances, the 
maximum levels of mixing are achieved for a value of AR Q /(b£/b o ) « 8. 

The results, as depicted in Fig. 4-18, can be formulated into a design criteria 
for rectangular unlike doublets for optimum mixing. These criteria are embodied 
in the following equations: 

Design for optimum mixing level 

<f> = 1.0 (4-7) 

and 

AR o /(b f /b o ) 3 = 8 (4-8) 

for w q = w^ only 

Of the two criteria, Eq. 4-7 is by far the most stringent. Deviations from <j> = 1 

can produce large reductions in the level of mixing. It can be seen from Fig. 

4-18, however, that for rounded entrance rectangular unlike doublets there is a 

great degree of latitude in the value of the constant assigned to Eq. 4-8. This 

constant can vary from 3 to 50 with no significant reduction in level of mixing. 

With sharp orifice entrances, the maximum obtainable mixing level is only a few 

3 

percent lower, but the sensitivity of off-optimum values of AR q / ( b^/b^ is sub- 
stantially greater. 
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L7. Rectangular Unlike- Doublet Mixing as a 
Function of b f /b and AR q for 4> = 1*00 
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Figure 4-18. Mixing Uniformity for Rectangular Unlike Doublets as a Function of 
AR 0 /(b f /b 0 ) 3 (Sharp and Rounded Entrance Data Shown) 




5.0 APPLICATION OF RESULTS TO INJECTOR 
DESIGN AND GENERAL CONCLUSIONS 


5.1 GAS/LIQUID ELEMENTS 

The results that have been presented for the gas/ liquid rectangular concentric 
tube injector elements suggest that the introduction of rectangular shape for 
injector design would be most beneficial for those applications in which high 
mixture ratio (MR > 4) and/or low gas momentum (p^,V^< 50 psi) are imposed.* For 
these applications, the rectangular shape can be employed to significantly im- 
prove mixing. However, care must be taken to ensure that the particular thrust 
chamber is not vaporization performance limited. (It has been shown that an in- 
creased aspect ratio produced an increased droplet diameter--see Fig. 3-43, An 
aspect ratio of at least 6:1 must be employed to achieve dropsizes that are com- 
parable to a circular element of equivalent injection areas.) 

In general, no significant additional mixing advantage is achieved for aspect 
ratios greater than 3:1 for zero center-post recess. However, higher aspect 
ratios do appear to offer added mixing quality with recessed center-posts (see 
Fig. 3-28). 

Aspect ratios greater than 6:1 are probably not practical from a fabrication 
standpoint. However, aspect ratios on the order of 6:1 are attractive as they 
provide the desirable mixing advantages and, at the same time, offer no apparent 
atomization penalty (see Fig. 3-43). 

The limited data obtained under this study suggest the guidelines listed in this 
section (see Table 5-1). 


*A rocket engine that operates with gaseous hydrogen at 800 psi will produce a 
value of P G V G 2 = 58 for V Q = 1000 ft/sec. 
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TABLE 5-1. SUGGESTED DESIGN GUIDELINES FOR APPLICATION OF 
RECTANGULAR CONCENTRIC TUBE INJECTOR ELEMENTS 


Quantity 

Target Valves 

Range of Application 

w /w >4.0 

Liquid Port Aspect Ratio 

l ~ 3:1 (mixing limited) 


X l^l ~ (mixing and vapori- 

1 zation limited) 

Center-Post Recess 

r /y l = 1.0 

Gas Port Aspect Ratio 

x 3 /y 3 =2.1 for x 1 /y 1 = 3.0 
Xj/y^ Unknown for x^yj “ 6 


5.2 LIQUID/LIQUID ELEMENTS 

Design guidelines and conclusions concerning optimum mixing for unlike triplets 
and doublets as well as some guidelines for triplet atomization design are pre- 
sented in this section. The mixing guidelines may be employed directly to in- 
jector design. However, caution must be exercised with the atomization criteria. 
The triplet atomization criteria obtained in this study are applicable to laminar 
flow only in orifices with rounded entrances. In addition, property corrections 
from wax to the desired propellant must be employed. It is interesting to note 
that for the fuel component (UDMH) , velocities up to 50 ft/sec (in a 0.070-inch- 
diameter orifice) or 50 ft/sec (in a 0.020-inch-diameter orifice) produce Reynolds 
numbers less than 10^. Thus, actual rocket applications will be concerned with 
laminar jets if the entrances to the orifices are rounded. 

5.2.1 Triplet Mixing Design Guidelines 

The design guidelines for triplet elements having circular orifices are quite 
simple. Basically, the element design should be configured such that the 
following criteria are met (see Table 5-2). 
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TABLE 5-2. TRIPLET INJECTOR DESIGN CRITERIA FOR MIXING 


For d 1 /d 2 > 1.0 

0 = 1.00 

For d x /d 2 <1.0 

0 = 0.714 


For a given propellant combination and specified mixture ratio, the value of 0 may 
be computed with the following relationship: 


Given 

Pj and p 2 

and mixture ratio = w^w^ 



[5-1) 


Since the value of 0 is specified for optimum mixing, the following equation 
represents the design criteria for triplets: 


d i s/i pTpi 

d 2 V2 0 MR 2 


[5-2) 


1 

Once the value of has been established (iteration may be required to fix the 
appropriate value 0 ? 0), the maximum level of mixing is determined (see Fig. 4-6C) . 


5.1.1. Unlike-Doublet Mixing Design Guidelines 


The design guidelines for optimum mixing are essentially identical for both cir- 
cular and noncircular orifice unlike doublets. That is: 

0 = 1 
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However, the noncircular orifices provide an additional degree of design flexi- 
bility. The importance of this additional flexibility is illustrated in the 
following discussion. 

The specification of a propellant combination and an operational mixture ratio, 
when combined with the criterion, 0 = 1, yields the following relationships: 

• for circular orifices 



• for rectangular orifices 


(5-4) 




(5-5) 


For the circular orifice, this is the end of the line. The mixing level associ- 
ated with the particular orifice diameter ratio must be accepted (see, for example. 
Ref. 3). This limitation is due to the fact that the aspect ratio for a circle is 
unity and does not provide an additional degree of freedom. 

However, for the rectangular elements, the value of the aspect ratio of the oxi- 
dizer remains a free variable, even after the ratio b f /b Q has been stipulated 
according to Eq. 5-5. The optimum value of this variable can be computed from 
the following relationship to achieve a design capable of the maximum levels of 
mixture ratio uniformity. 


(AR ) . 

o' opt 




(P c /P f ) 


MR 


(5-6) 
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An approximate dropsize relationship has been presented for triplet elements 
with circular (rounded entrance) orifices: 


D = 2.15 x 10 S V“- 77 



(5-7) 


where 

V. ' = ft/sec 
J 

Dj = inches 
D = microns 

This relationship may be employed to compute the approximate mass median wax 
droplet diameter produced by either the central orifice or the outer orifices 
for laminar flow. However, there is a strong property correction that must be 
applied to the wax droplet diameter. This correction has not been established 
for laminar jet flow with triplet elements. Yet, some appreciation for the 
relative importance of the changes of the injector design variable can be real- 
ized. The relationship for droplet diameter can be employed to test the sensi- 
tivity of dropsize to injection velocity (Ap) , orifice length, and jet diameter. 

For example, the relationship suggests that V. and D. are of equal import while 

L J J 

the orifice length (^) is less important (or influential) . 
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6.0 RECOMMENDATIONS FOR FUTURE WORK 


6.1 GAS/LIQUID WORK 

In the areas of gas/liquid rectangular and circular concentric tube mixing and 
atomization, further cold- flow experimentation is required with variation of area 
ratio. This will aid in generalizing the results that are presently restricted 
according to Eq. 3-8. Some work in this area has been reported in Ref. 7. However, 
these results were limited to a specific application. More general parametric data 
are needed. 

Further data concerning the effect of gas port aspect ratio, x^/y^, upon the mixing 
characteristics for rectangular concentric tube elements (RCTE) , are required. 

The results presented in this report show this to be a potentially valuable param- 
eter for the optimization of mixing, and also for control of the spatial mixture 
ratio distribution produced by injector elements. The latter feature could have 
valuable application with regard to injector/thrust chamber compatibility. 

The work with RCTE’s should be carried to its natural conclusion- -the design and 
hot-fire demonstrations of a large-scale injector designed for high mixture ratio 
(e.g., LOX/GH^ at MR = 6) and low gas injection momentum. This type of study 
could also include a variation of the orientation of the peripheral elements to 
evaluate the advantage of RCTE's for improved injector/ chamber compatibility, and 
perhaps stability. 

6.2 UQUI D/LIQUID WORK 

An evaluation of the mixing characteristics of rectangular unlike doublets with 
unequal facing widths (i.e., w q f w £ , see Fig. 4-12) is strongly recommended. The 
additional flexibility of the parameter w q /w £ , when added to Eq. 4-7 and 4-8 could, 
perhaps, greatly increase the maximum levels of single- element mixing available 
with uni ike- doublet elements. This is an extremely valuable goal due to the fact 
that the simplicity of the unlike-doublet element makes it an attractive candidate 
for injector design. 
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In conjunction with the mixing evaluation, a detailed study of reactive stream 
separation ("blowapart") should be conducted. If the unlike doublet is to be 
employed with hypergolic propellants, this problem must be address and definitive 
design guidelines for the avoidance of blowapart established. 
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APPENDIX A 


GAS/LIQUID PRESSURIZED COLD-FLOW MIXING FACILITY* 
MEASUREMENT SYSTEM 


To characterize the spray fields generated by gas/liquid rocket motor injectors, 
a system was developed to determine local values of gas and liquid mass flux 
under pressurized conditions. Knowing local values of gas and liquid mass lux, 
the "mixedness" (i.e., E m ) of the two-phase spray field was determined. 

A schematic of the complete measurement system for the characterization of dense 
gas/liquid spray fields is shown in Fig. A-l. Mixing experiments were performed 
with both single-element and multi-element gas/liquid injectors. The apparatus 
consisted essentially of a pressurized test section in which the deceleration 
probe (described in subsequent paragraphs) was positioned in r-0 coordinates. 

The system contained several "water traps" to ensure that water, which can accumu- 
late during extended test periods, does not plug critical pressure lines in the 
system. A photograph of the system, which is located at the Combustion and Heat 
Transfer Laboratory, is shown in Fig. A-2. 

A problem associated with the characterization of spray fields generated with 
single-element injectors was the suppression of the flowfield recirculation 
caused by the injection of high-velocity streams into a finite closed volume. 

High levels of recirculation within the test section precluded the accurate de- 
termination of the gas-phase flowfield. To eliminate flowfield recirculation, a 
low-velocity uniform "base bleed" flow surrounded the single-element injector. 

The values of base bleed velocity are in accord with the CrayaCurtet criterion 
for the elimination of recirculation (Ref. A-l). 

One additional problem associated with the characterization of spray fields 
generated with single-element injectors was the determination of the local mass 
flux of the injectant gas. As the gas/liquid flow field moves through the 

*Appendix A taken directly from Ref. A-2. 
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Figure A-2. Pressurized Gas/Liquid Mixing Facility 



surrounding environment on its way from the injector to the probe tip, much of 
the gas in the environment (base bleed) was ingested into the flowfield. There- 
fore, the gas flowrate measured by the probe was composed of gas that was injected 
and gas that was ingested. These two components were segregated to determine the 
injected portion (for calculation of E m ) . To accomplish this, nitrogen and oxygen 
were mixed in known quantities and injected through the element while pure nitrogen 
was used as a base bleed fluid. The sample that arrived at the probe tip was 
processed through a Beckman oxygen analyzer (Fig. A-l). The concentration of the 
oxygen in the sample was then used to determine the concentration of the injected 
gas phase. 

TWO-PHASE DECELERATION PROBE 


The two-phase deceleration probe, which was used for the determination of local 
values of gas and liquid mass flux, was developed under Contract NAS3-12051 (see 
Ref. A-2). Only a brief description of the probe will be presented herein; for 
a more detailed description see Ref. A-3. 

The concept of a deceleration probe for the measurement of gas-phase stagnation 
pressures, for the determination of gas velocity,, in two-phase flowfields was 
first introduced by Dussourd and Shapiro (Ref. A-4) . However, the referenced 
probe design was operated only in low mass flux ratio (particle flowrate/gas 
f lowrate < 0 . 2) flowfields. In addition, the probe design of Dussourd and 
Shapiro was utilized only in ambient pressure flowfields. 

The deceleration probe which was developed at Rocketdyne was utilized both for 
the determination of gas-phase stagnation pressures and local liquid mass flux. 
The probe has been demonstrated to operate successfully in high mass flux ratio 
(liquid mass flux/gas mass flux 0.2 to 20) two-phase flowfields (Ref. A-5) . In 
addition, measurements have been made in dense gas/liquid flowfields at static 
pressures up to 500 psia. 


V/ 
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A schematic of the deceleration probe, termed the "concentric tube two-phase 
impact probe", is presented in Fig. A-3. The probe was constructed of two con- 
centric tubes (A and B) with a specially designed tip attached to tube B. The 
tip was designed to prevent the passage of water (termed flooding) into the 
annulus formed by tubes A and B when the probe is utilized in high mass flux 
ratio flowfields. The problem of flooding is a serious limitation of the probe 
design described in Ref. A-4. 

The operating principle for the determination of the gas-phase stagnation pressure 
by the concentric tube two-phase impact probe is illustrated in Fig. A-4. Basi- 
cally, the intent is to decelerate the gas and measure the gas-phase stagnation 
pressure in a manner that minimized momentum exchange from the condensed phase 
upstream of the measurement location. Particles. (or droplets) and gas (each at 
their own velocity) encounter the probe tip. The gas phase is stagnated at the 
probe tip where the pressure is approximately equal to the gas-phase stagnation 
pressure. Deviation from true gas-phase stagnation pressure is due to momentum 
exchange between the particles and the gas in the near flowfield of the probe 
tip (termed overpressure error). A particle, due to its higher inertia, passes 
through the probe tip and is decelerated to zero velocity in the stagnation 
chamber formed by tube A. However, due to momentum exchange between the particles 
and the stagnated gas, the particles decelerate in the probe tip to some extent 
over the distance X (see Fig. A-4). The gas-phase stagnation pressure Po,gas» as 
measured in the probe annulus, is greater than the gas-phase stagnation pressure, 
p o,tip- The difference between the two aforementioned pressures can be made small 
if the distance X is minimized. However, the total overpressure error (due to 
particle/gas momentum exchange both near and within the probe tip) can be deter- 
mined by proper calibration of the probe in known two-phase flowfields (Ref. A-3). 

The gas-phase stagnation pressure was measured in the 360-degree annulus at the 
probe tip rather than at a single point as was done on the probe described in 
Ref. A-4. Tests conducted with both the subject probe and one of the Ref. A-4 
designs in high mass flux ratio (>3) flowfields demonstrated that the concentric 
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Figure A-4. Principle o£ Operation 


tube probe eliminated the pressure oscillations encountered with a probe of 
Ref. A-4 design. In addition, the concentric tube probe design avoids the need 
for an external pressure line at the probe tip that could create flowfield dis- 
turbances. The gas-phase stagnation pressure in the probe annulus was measured 
with an MKS Baratron Type 77 electronic pressure meter. 

Local values of liquid mass flux were determined by capturing a liquid sample 
in the probe stagnation chamber (Fig. A-2) for a known time interval. However, 
since the droplet capture efficiency of the probe design is less than 1 (aiO.95) 
small corrections must be made to the captured liquid mass to determine a value 
of local liquid mass flux. Values of the probe capture efficiency were deter- 
mined by calibration of the probe in known two-phase flowfields. 

A photograph of the two-phase probe and its traversing mechanism is shown in Fig. 
A- 5. The probe is shown located in the rotatable portion of the test system. 
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APPENDIX B 


GAS/LIQUID PRESSURIZED COLD-FLOW ATOMIZATION FACILITY 


TECHNIQUE 

To simulate the dynamic injection conditions of hot-fire propellant systems, the 
atomization experiments were conducted at simulated gas-phase densities. In this 
case, gaseous nitrogen and molten wax (Shell-270) were used as nonreactive fuel/ 
oxidizer simulants, respectively. In this molten wax technique, the wax droplets 
freeze before collection and are subsequently subjected to sieve analysis after 
drying. 

A series of 23 standard sieves ranging from 53 to 2380 microns (53 to 2380 x 10 
meters) was available for sample sieving. For any particular sample, only 12 of 
the sieves were used, the particular sieve sizes that were used depended on the 
anticipated size range. The particular wax sample (~10 grams) was shaken on a 
RO-TAP automatic shaker for 30 minutes. After the sieving operation was completed, 
the mass of particles on each sieve was weighed on an electric balance. A total 
recovery of 96 to 98 percent of the mass originally introduced into the sieves 
was possible. The mass fraction of sample on each sieve was then plotted as a 
function of sieve size to determine a mass median drop size. 


EXPERIMENTAL FACILITY 

A schematic of the pressurized atomization facility that was employed in these 
studies is presented in Fig. B-l. The system consisted of a 600-gallon cylincrical 
tank in which a single-element injector model was mounted. Molten wax and heated 
GN were supplied to the injector from a system that was heated with circulating 
hot oil. All lines and valves in the wax supply were oil jacketed to prevent wax 
freezing. Heated GN^ was supplied to the injector at a temperature above the melt- 
ing point of the wax («200 F) so that the resulting wax droplets did not freeze 
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Figure B-l. Schematic of Pressurized Atomization Facility 



before the completion of the liquid jet breakup and atomization processes. A water 
flush was supplied to the inner tank walls during testing to ensure that molten wax 
did not adhere to the tank wails. A photograph of the pressurized atomization 
facility is shown in Fig. B-2. 
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Figure B-2. Pressurized Atomization Facility 
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APPENDIX C 


LIQUID/LIQUID MIXING FACILITY 


COLD-FLOW FACILITY 


Flow System 


The basic components of the flow system are shown in Fig. C-l. The system contains 
two high-pressure (1000-psia maximum rated pressure) supply tanks. Each are re- 
motely pressurized. The propellant lines are stainless-steel tubing. Pneumatic 
(Annin) valves are used for tank shutoff and main valves. 

Flow system instrumentation consists of four Taber "Teledyne" series-bonded strain 
gage pressure transducers and two Fischer -Porter turbine-type volumetric flow- 
meters. Measurements are made of both propellant tank pressures and propellant 
injection pressures. These measurements are recorded on Dynalog direct-inking 
graphical recorders. The volumetric flowmeter signals are recorded on a CEC 
multichannel oscillograph. 

Cold-flow propellant simulants are trichloroethylene and water, which simulate 
the oxidizer and fuel, respectively. These simulants were chosen on the basis of: 
(1) availability, being employed on related programs at Rocketdyne using the same 
facility, (2) ease of handling, and (3) maintenance of the oxidizer/fuel 
immiscibility . 

Collection System 

The specific details of the collector are illustrated in Fig. C-2. As can be 
noted, the tubing slants outward from the collection plane to a 7- by 7-foot 
base. The base is 1/2-inch aluminum plate and separates the upper portion of 
the assembly from the Pyrex tube racks. Beneath the aluminum plate is a cart that 
houses the tube racks. The cart is mounted on wheels so that the entire tube 
matrix is easily removed from under the collector and rolled to the measurement 
station. 
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Figure C-l. Schematic Drawing of the Flow System for the Cold-Flow 
Mass and Mixture Ratio Distribution Study 
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Figure C-2. Injector Spray Collection System 
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EXPERIMENTAL TEST PROCEDURES 


The procedure for conducting a cold-flow test are as follows. The fuel and oxi- 
dizer simulant tanks are pressurized to give the desired flowrates. The main 
valves are opened and, after the injection pressures become steady, the shutter is 
opened for a selected time interval and then closed. The main valves are then 
closed to conclude the test run. All tests are conducted with the elements cen- 
tered above the collector at a distance of approximately 3 inches. This distance 
was chosen since previous analytical and experimental data indicate that it repre- 
sents a good approximation of the primary propellant mixing region during combustion. 

The data recorded for each test included oxidizer and fuel simulant flowrates, 
injection pressures, flow duration, and the volume of oxidizer and fuel simulants 
in each of the collection tubes. 

The individual volumes are determined by a volumetrically calibrated metal strip 
or graduated cylinders. The metal strip is a thin metal scale with a scribe 
mark at 10-milliliter increments. This strip is inserted into the test tube and 
the volume of oxidizer simulant and the volume of fuel simulant read directly. 

For tubes containing insufficient liquid quantities for accurate measurement, the 
volumetric measurements are obtained by use of graduated cylinders. 

DATA ANALYSIS 

The collector matrix data and the other recorded data are processed by computer 
to produce the following output: mass of oxidizer simulant, mass of fuel simulant, 

mixture ratio, and mass fraction for each tube. The mixing factor (E ) , cold flow 

m , 

c* efficiency * center of collected mass (row and column), and percentage 

of the injected mass collected are computed. 
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APPENDIX D 


LIQUID/LIQUID ATOMIZATION FACILITY 


COLD- FLOW FACILITY 

The wax flow facility used for the liquid/liquid dropsize measurements is shown 
in Fig. D-l and D-2. The overall system consists of wax and water supply tanks 
immersed in hot oil bath container (same system employed for gas/liquid studies) 
and a particle collector that catches the frozen wax particles. Instrumentation 
requirements of pressure, flowrate, and temperature are provided by strain gage 
transducers, turbine flowmeters, and iron constantan thermocouples, respectively. 
Each wax and water tank has an independent pressurizing and vent system. Also, 
each product outline has three flowmeters, thermocouples, and hand shutoff valves 
arranged in parallel so a wide range of flowrates can be obtained. The hot oil 
bath (shown schematically in Fig. D-2) is heated by means of a 30-kw thermostat- 
ically controlled heater. An electrically operated pump circulates the oil from 
the oil bath container through the heater and back again to ensure uniform tem- 
perature. Also, hot oil is forced through jacketed run lines and valves to ensure 
that the wax does not freeze in the feed lines. 

The particle collector (Fig. D-l) is an 18- by 50-foot, epoxy-coated, wooden plat- 
form that is located within a building. This shields the collection area from 
wind currents that might cause the smaller particles to be blown away. When the 
impact surface is washed down, the slope of the collector causes the wash water 
to be directed into a particle catch basin. The catch basin has several baffles 
to ensure that none of the wax particles are washed overboard. 

Experimental Procedure 

The experimental procedure for droplet size measurement is as follows: 

1. The proper injector configuration is installed on the wax facility in 
such an orientation that the wax spray, after freezing during its bal- 
listic trajector, strikes the particle collector. 
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Figure D-l. Open-Air Atomization Facility 



Figure D-2. Schematic of Hot-Wax Flow Facility 
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2. The electrical oil heater and pump are turned on to bring the propellant 
simulant tanks and run lines up to 210 F. 

3. After all parts of the system are thoroughly heated and instrumentation 
requirements checked, the run tanks are pressurized. 

4. With the piston -operated shutter in the up position, the test is initiated 
by actuating the main pneumatic shutoff valves. When the flowrates and 
injection pressures reach a steady condition, the shutter is actuated 

and the wax particles are allowed to spray onto the particle collector. 

The use of the shutter minimizes the influence of start and stop tran- 
sients on the size distribution of the collected particles. 

5. The injector flow is continued for approximately 10 seconds. The shutter 
is then actuated to the up position and main shutoff valves closed. 

6. The tanks are then vented and systems secured. 

7. The particles are then washed from the collector into the catch basin, 
where they are scooped from the surface of the water and placed in a 
plastic bag for temporary storage. 
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APPENDIX E 


HOT-FIRE EXPERIMENTAL FACILITY 


The gas/liquid hot-fire experimentation was conducted during Phase IV of this 
program on test stand Lima at the Propulsion Research Area (PRA) located at the 
Rocketdyne Santa Susana Field Laboratory Facility. A schematic of that test 
stand is presented in Fig. E-l with the critical instrumentation keyed, numer- 
ically, to the instrumentation list presented in Table E-l. 

The recording method for each parameter is noted in Table E-l. All of the param- 
eters critical to performance calculation were recorded by the Beckman automatic 
data acquisition system. 

Pressures were measured with Taber strain gage type pressure transducers. All 
temperatures were measured with iron- constant an thermocouples with output in 
millivolts. The liquid oxygen and the coolant water flowrates were measured with 
Fischer-Porter turbine flowmeters. The hydrogen gas flowrate was metered with a 
Flowdyne sonic venturi. 

As shown in Fig. E-l, there were four primary systems involved in the operation 
of test stand Lima. There were the primary propellant systems, liquid oxygen and 
gaseous hydrogen, a water system for thrust chamber coolant, and a gaseous fluorine 
system for ignition. (GF^ was contained in a single, 400-psi maximum commercial 
gas cylinder.) 

All valve opening and closing was controlled by a programmed, automatic timing 
system. The start and shutdown sequence is listed in Table E-2, 
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Figure E-l. Test Stand Lima Schematic, 1/2 lb/sec 
LOX/GH 2 Facility With GF 2 Ignition 
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TABLE E-l. LIMA STAND INSTRUMENTATION LIST 




Type or 
Range 


Recording 

7 

> 


Parameter 

Beck. 

Osc. 

DIGR 

1 . 

Pressure, LOX Tank 

0-3000 psi 



X 

2. 

Temperature, LOX Flowmeters 

I.C. 

X 



3. 

Flowmeter 1-LOX (5806A-3173-A2) 

Turbine 

X 

X 

X 

4. 

Flowmeter 2-LOX (5808A-3673-A2) 

Turbine 

X 

X 

X 

5. 

Pressure, Upstream LOX Venturi 

0-2000 psi 

X 


X 

6. 

Temperature, Upstream LOX Venturi 

I.C. 

X 



7. 

Pressure, LOX Injection 

0-2000 psi 

X 

X 

X 

8. 

Temperature, LOX Injection 

I.C. 

X 


X 

9. 

Pressure 1, Upstream GH 2 Venturi 

0-2000 psi 

X 


X 

10. 

Pressure 2, Upstream GH^ Venturi 

0-2000 psi 

X 



11. 

Temperature, Upstream GH 2 Venturi 

I.C. 

X 



12. 

Pressure, GH 2 Injection 

0-2000 psi 

X 

X 

X 

13. 

Pressure, H 2 0 Tank 

0-3000 psi 



X 

14. 

Temperature, H 2 0 Heater Line 

I.C. 



X 

15. 

Flowmeter 1, H 2 0 (Total) 

Turbine 

X 

X 

X 

16. 

Pressure, H 2 0 Inlet 

0-2000 psi 

X 


X 

17. 

Delta Temperature, Coolant H 2 0 

3-Leg 

Thermopile 

X 


X 

18. 

Flowmeter 2, H 2 0 (Chamber) 

Turbine 

X 

X 

X 

19. 

Pressure, Chamber H 2 0 Outlet 

0-2000 psi 

X 



20. 

Temperature, H 2 0 Inlet 

I.C. 

X 


X 

21. 

P -1 (Injector Face) 

0-1000 psi 

X 

X 

X 

22. 

P -2 
c 

0.1000 psi 

X 

X 

X 

23. 

P -3 
c 

0-1000 psi 

X 

X 

X 

24. 

Pressure, GF 2 Supply Line 

0-1000 psi 



X 

25. 

Temperature , LOX Line 

I.C. 

X 



26. 

Pressure, LOX Flowmeters 

0-2000 psi 

X 
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TABLE E-2 . OPERATIONAL SEQUENCE FOR LIMA STAND 


• Pressurize H 2 0 and LOX tank and pressurize GH 2 

• GN^ purges on 

• Water coolant main valve open 

• GH 2 main open 

• GF 2 main open 

• LOX main open 

• GF 2 main closed 

• LOX main closed 

• GH 2 main closed 

• Coolant main closed 

• GN 2 purges off 

• Vent all tanks and loaders 


loader 


R-9270 

E-4 



APPENDIX F 


COMBUSTION MODELS 


VAPORIZATION-LIMITED COMBUSTION 

The vaporization- limited combustion model formulation is based on the development 
of mathematical expressions for the various physical processes involved in the 
combustion of liquid droplet sprays in a bipropellant liquid rocket engine. The 
model considers propellants to be injected as sprays containing ranges of discrete 
droplet size groups, each possessing a given average diameter. The total spray 
mass is distributed among the various groups according to an experimental mass 
distribution function. 

Of central importance in the model is the solution of the individual droplet burn- 
ing rates, which are assumed to be limited by diffusion. Analysis of the dynamic 
behavior of single droplets is justified on the basis that the volumetric flowrate 
of liquid propellants into the downstream region is only about 1 to 2 percent of 
that of the combustion gases and, therefore, that the likelihood of droplet colli- 
sions or interference with one another is negligibly small. Under rocket condi- 
tions, in the uniform mixing zone, droplets are spaced on the order of 2 to 3 
diameters apart, while the vapor film thickness is on the order of 5 to 15 percent 
of the droplet diameters. As a result, each droplet is considered to be immersed 
in an infinite combustion gas medium. 

The calculation of single droplet evaporation is based on a spherically symmetric 
model of simultaneous heat transfer to, and mass transfer from a liquid sphere. 

The liquid droplet temperature is assumed to vary with time, but to be uniform 
through the drop. Forced convection and resultant nonspherical transfer processes 
are accounted for through empirical Nusselt number correlations for both heat and 
mass transfer. 
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In evaluating the convective contribution, relative gas-to-droplet velocity is 
required. Droplet velocities are obtained from a drag relationship for evaporating 
spheres. A composite form of the drag coefficient for accelerating spheres that 
account for droplet flattening is employed. 

Compressible gas dynamics are accounted for with area changes corresponding to 
chamber geometry. The droplets are treated as point sources (or sinks) of fuel 
(mass), oxidizer (mass), momentum, and energy with local transport rates obtained 
by summing the contributions of all droplets at any given location in the chamber. 

The gas-phase energy equation is simplified normally by the assumption that the 
composition and stagnation temperature are the equilibrium values for the gas- 
phase oxidizer, fuel mixture ratio, and the chamber pressure. Other gas proper- 
ties (static temperature, density, etc.) are evaluated from the respective stagna- 
tion values by applying the local Mach number to the frozen isentropic expansion 
equations . 

The model is solved in numerical form by high-speed digital computers. It re- 
quires input of the "upstream boundary condition,” which completely describes the 
initial conditions of spray (dropsize distribution, drop velocities, and tem- 
perature) and gas (composition, flowrate, and pressure) at the location where com- 
putation is started. Chamber geometry must also be specified. 

Solution proceeds in a stepwise manner moving downstream to the nozzle throat. 

At each step, interphase transport of mass, momentum, and energy is evaluated from 
the transport equations previously described with subsequent solution of gas-phase 
equation of state and continuity, momentum, and energy balances. This results in 
a description of droplet group diameters, velocities, and temperatures as well as 
gas composition, velocity, and pressure at the new location. 

This "marching technique" proceeds into the nozzle up to the geometric throat, 
where it is necessary to satisfy the downstream boundary condition of sonic gas 
velocity. If the throat Mach number deviates from unity by more than a pre- 
selected tolerance, iteration is required whereby propellant flowrates are adjusted 
and the entire calculation repeated. In practice, convergence of this iteration 
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is rapid and a solution is readily obtained. The general validity of the analy- 
tical results is determined to a major extent by the accuracy of the input spray 
description. Vaporization rate limited c* efficiency is computed from the follow- 
ing equation: 


,vap 


R] ft] 


(P-X) 


where 

Wg = flowrate of burned gas at the geometric throat 
Wj = injection flowrate of fuel plus oxidizer 

c* B = theoretical c* corresponding to the composition of the burned gas 
at the geometric throat 

c*j = theoretical c* corresponding to the injection mixture ratio of 
liquid fuel and oxidizer 


MIXING LIMITED COMBUSTION 


Over the past 15 years, mass and mixture ratio distribution uniformity ("mixing") 
has been extensively studied both analytically and experimentally. Experimental/ 
analytical correlations demonstrate quantitatively that high combustion efficiency 
in rocket engine thrust chambers occurs only when the initial local mixture ratio 
distribution is at, or near, the target chamber mixture ratio. This implies that 
the injector should provide a spray field having a uniform mixture ratio over the 
entire flow cross section. 

The sketch shown on the following page illustrates a typical curve of theoretical 
equilibrium c* as a function of propellant mixture ratio (oxidizer/fuel). Nor- 
mally the design operating point of overall injected mixture ratio falls close to 
the peak, and any maldistribution of propellant mixture ratio results in a loss 
in overall c*. An analytical model has been developed at Rocketdyne to relate 
these maldistributions to an attendant loss in c* efficiency. The development of 
this model is outlined in the following paragraphs. 
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Wrubel (Ref. F-l) describes an analysis of mixing losses whereby the flow is hypo- 
thetically subdivided into "i" stream tubes, each containing propellant at some 
mixture ratio, which is uniform within that stream tube. No mass or energy is 
considered to cross stream boundaries. Propellant vaporization, mixing, and com- 
bustion are treated as being complete upstream of the start of nozzle convergence. 
Within the nozzle the flow is handled as being one-dimensional and isentropic. 

At each axial station the static pressure is considered uniform for all stream 
tubes. In addition, boundary layer effects are neglected. The resulting equa- 
tion relating the mixing limited c* efficiency to the local mass and mixture ratio 
distribution is: 


n • 

mix 





c* 

theo 



(F-2) 


Here A^ ./A*, is the ratio of the cross-sectional area of the ith stream tube at 
t,i 1 

the minimum chamber area to its area at the point it becomes sonic. For most 
cases of interest the specific heat ratios y^, are all of similar value so that 
shifts in location of the sonic condition from the geometric throat will be small 
and the preceding equation is closely approximated by: 


n . = E 

mix 

1 


w. 

1 


w c* 

T theo 


(F-3) 


R-9270 

F-4 



where the effective c* is simply a weighted average of the local c* for the indi- 
vidual stream tubes. For any given propellant mixture ratio distribution, Eq. F-3 
provides a simple means of determining c* efficiency loss due to "mixing." 

Most investigators agree that distributions developed by spray mixing near the 
injector will not be appreciably changed downstream by turbulent mixing of the 
gases. As a consequence, if the initial spray distribution formed by an injector 
can be experimentally determined, (h c *) m ^ x can be computed by using Eq. F-3. 

Rupe (Ref. F-2) introduced a term, commonly known as E m , which is an index of mixing 
uniformity: 



R - ratio of total oxidizer mass to total oxidizer and fuel mass 

r i = ratio of oxidizer mass to total oxidizer and fuel mass in an 

individual stream tube for r^ < R 

r. = ratio of oxidizer mass to total oxidizer and fuel mass in an 

individual stream tube for r. > R 

1 

The factor E^, is not uniquely defined by C 7 ! corres P on ^ enc 6 is strongly 
affected by the propellant combination and the nominal mixture ratio. This term 
is employed to describe the average mixing uniformity of a given spray field. 
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APPENDIX G 


PROPERTIES OF FLUIDS EMPLOYED IN PHASES III AND IV 
COLD- FLOW AND HOT-FIRE PROGRAMS 

The properties for the various fluids listed in Table G-l are presented in 
Fig. G-l through G-3. 


TABLE G-I. KEY TO PROPERTY CHARTS 
(Numbers refer to figure numbers) 



Property 

Viscosity 

1 

Surface Tension 

Acoustic 

Velocity 

Liquid Oxygen 


G-2 


Shell Wax (270) 

1 ;:' 

G-2 


Water 

G-l 

G-2 


Gaseous Nitrogen 

G-l 


G-3 

Gaseous Hydrogen 

G-l 


G-3 
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Figure G-l. Viscosity of Liquid Oxygen, Wax (Shell 270), Water, 
Gaseous Nitrogen, and Gaseous Hydrogen as a Function 
of Temperature 
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Figure G-3. Acoustic Velocity of Gaseous Nitrogen 
and Gaseous Hydrogen as a Function 
of Temperature (Low Pressure) 
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APPENDIX H 


DISCUSSION OF THE CENTERLINE MOMENTUM RATIO PARAMETER, 0 


INTRODUCTION 

Mixing characteristics of two impinging jets of dissimilar fluids were first 
studied by Rupe (Ref. H-l) in 1953. Rupe evaluated the effect of density, veloc- 
ity, diameter, and impingement angle on the mixing uniformity and found the re- 
sultant mixing characteristics strongly dependent upon these variables. In 
particular, Rupe has shown that a given unlike impinging doublet injector (cir- 
cular orifice) produces optimum mixing uniformity when the product of density, 
velocity-squared and the diameter for each of the jets are equal, and that this 
optimum level is a function of geometry. It is not clear from Rupe's work whether 

the variables p, U, D, are all independent, composed of a dynamic pressure term 

2 2 
(pU ) and a characteristic dimension, or a single variable(pU D) . Over the past 

15 years Rupe's result has been routinely applied to design injectors without a 
clear understanding of the significance of the fluid dynamics defined by his 
relationship. From an applications standpoint this is understandable and pre- 
sented no design problems for circular orifices. However, reformulation for non- 
circular designs results in several possible forms of Rupe's criteria, depending 
upon the specification or grouping of the variables. Consequently, application 
of Rupe's criteria to noncircular orifice designs requires understanding of the 
mechanisms controlling mixing. 

DISCUSSION 

Based on Rupe’s work, a dimensionless quantity which his experiments have shown 
important to mixing is : 

P 1 U l 2 °1 /P 2 U 2 2 °2 (H_1) 
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Inspection of the physical significance of the variables contained in Eq. H-l for 
circular geometry shows that the variables can be reformulated in terms of 
dynamic pressure (0) 

0(D f /D Q ) =1.0 (H-2) 

or momentum ratio of the jets (M^/M q ) , and a diameter ratio: 

V M 1 (VV = 1,0 (H ' 3) 

It becomes obvious that neither of the above physical interpretations result in a 
meaningful dimensionless group since they still result in two physical dimension- 
less ratios. 

Since for circular orifices the above equations are dependent, application of 
either definition (Eq. H-5 or H-6) will lead to the proper design values. However, 
if the criteria are extended to other orifice geometries wherein the equations 
become independent (i.e., noncircular orifices), then the physical significance 
of Eq. H-l must be known to select the proper form. As an example, a set of rec- 
tangular orifices impinging at 60 degrees, having equal widths, result in the 
following equations, which are equivalent to Eq. H-2 and H-3: 


i 

W 

— T 

b i I — H b 2 I— 


0 <V D oVdrauli= ' CP 1 U l 2/P 2 U 2 < Whydraulic 


(H-4) 
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(W'Whydraulic 


- <Pl u l 2 V p 2 U 2 2 VC57) 


(H-4) 


hydraulic 


In this example, the hydraulic diameter ratio has been substituted for the geo- 
metric term. Note that these equations are independent and, therefore, should 
result in differing correlations of data. In addition, since the physical signi- 
ficance of the original equations is not clear, it is possible that neither of 
these correlations will result in a meaningful description of the mixing process. 

Equation H-l can also be formulated in terms of the centerline momentum of the 
jets. It is easily shown that on the basis of centerline momentum, Eq. H-l for 
circular orifices is: 


_L 

°z 

T 





(a) CIRCULAR JETS 



dx-H 


2 

T 


(b) RECTANGULAR JETS 


0 = p f U/ D f dx/P o U 0 ‘ D o dx - p f U f D f /p 0 U 0 D o 


and for rectangular orifices is: 


0 = P, U 


. 2 V p 2 u 2 2 b ; 


CH-6) 


(H-7) 


Physical arguments, detailed in Ref. H-2, were employed to show that Eq. H-7 is 
the proper form for the parameter 0. 
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